IPNO/TH 95-45 
PURD-TH-95-05 



in 

On 



3 



X 



The Low Energy tttt Amplitude 
To One and Two Loops 



M. Knecht, B. Moussallam, J. Stern 

Division de Physique Theorique f\, Institut de Physique Nucleaire 
F-91406 Orsay Cedex, France 



>• ! and 
ON 

£2 '■ N. H. Fuchs 

O 

in 

; Department of Physics, Purdue University 

West Lafayette, IN 47907, USA 



Abstract 

The low-energy tttt amplitude is computed explicitly to two-loop accuracy in the chiral 
expansion. It depends only on six independent (combinations of) low-energy constants 
which are not fixed by chiral symmetry. Four of these constants are determined via 
sum rules which are evaluated using tttt scattering data at higher energies. Dependence 
of the low-energy phase shifts and of the threshold parameters on the remaining two 
constants (called a and f3) are discussed and compared to the existing data from K14 
experiments. Using generalised xPT, the constants a and j3 are related to fundamental 
QCD parameters such as the quark condensate (0|gg|0) and the quark mass ratio m s /fh. 
It is shown that forthcoming accurate low-energy tttt data can be used to provide, for the 
first time, experimental evidence in favour of or against the existence of a large quark- 
antiquark condensate in the QCD vacuum. 



* Unite de Recherche des Universites Paris XI ct Paris VI associee au CNRS. 



1 Introduction 



Current understanding of the mechanism of spontaneous breakdown of chiral symmetry 
(SB%S) in QCD is still lacking clear experimental support. The basic fact that in the limit 
m u = = m s = the chiral symmetry of the lagrangian is spontaneously broken down to 
S77(3)v is known to be a mathematical consequence of anomalous Ward identities quark 
confinement and the vector-like character of the theory ||. It just means that there are eight 
massless Goldstone bosons coupled to eight conserved axial-vector currents. In the real world 
where m u , rrid, m s are nonzero but much smaller than the typical mass scale Ah ~ 1 GeV of 
the first massive hadrons, the theoretical fact of SB^S and the Goldstone character of the eight 
lightest pseudoscalars are experimentally well founded. The chiral world, i.e. the first step in 
the expansion in powers of m quar k/A#, indeed strongly resembles the observed world. 

That much is known and is sufficient to establish a low-energy effective theory in which 
the role of effective degrees of freedom is played by weakly interacting Goldstone bosons ||. 
Equivalence of this effective theory with QCD is ensured step by step in a systematic expansion 
in powers of external momenta and quark masses, referred to as chiral perturbation theory 
[fj. ||] (%PT). The standard wisdom, however, involves much more than this established general 
framework of SB^S and of the corresponding low-energy effective theory since, in addition, it 
assumes a particular mechanism of symmetry breaking. What is commonly assumed is that in 
QCD, as in the Nambu-Jona-Lasinio model |§, SB^S is triggered by strong condensation of 
quark- ant iquark pairs in the vacuum. The quark condensate parameter 

Bo = ~ hm (1-1) 

is then assumed to play a dominant role in the description of symmetry-breaking effects 
induced by quark masses. The assumption of a strong quark condensation and its consequences 
have not so far been tested experimentally. Moreover, in a vector-like theory like QCD, a 
natural alternative scenario is conceivable. Goldstone bosons with a non-vanishing coupling 
F = lim m<juark ^o F-k to conserved axial currents can be formed even if quark condensation is 
marginal or even absent. For all these reasons an experimental probe of the mechanism of SB^S 
in QCD becomes of fundamental importance. It has been pointed out that the most promising 
framework for testing the strength of quark condensation is low-energy tttt scattering || ||. 
In the present work the theoretical basis for such a test is elaborated in detail, in view of two 
experimental projects to improve considerably the accuracy of existing low-energy tttc scattering 
data flO], fy]. The tttt scattering amplitude is worked out explicitly to two- loop accuracy in the 
general low-energy expansion, independently of any prejudice concerning the size of the quark 
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condensate. 

Before describing the content of this article in detail, let us briefly comment on the theo- 
retical possibility that in QCD, SB^S need not necessarily imply a strong quark condensation. 



This is illustrated by an analogy with antiferromagnetic spin systems ||12|| . In this case, the 
spontaneous breakdown of rotation symmetry does not necessarily yield a large spontaneous 
magnetisation (as it would happen for a ferromagnet). The latter can even vanish to the extent 
that the structure of the ground state reaches the Neel-type magnetic order. From renormalised 
perturbation theory, we are used to consider as "unnatural" a situation where a quantity would 
vanish (or be kept unusually small) if it is not forced to do so for reasons of symmetry. The 
above example of antiferromagnetism suggests that this reasoning need not necessarily apply 
to quantities of genuine nonperturbative origin such as spontaneous magnetisation or the quark 
condensate. Indeed, in QCD (as in any local relativistic field theory) the necessary and suf- 
ficient criterion for SB%S is the non-vanishing of the coupling Fq of the Goldstone bosons to 
the Noether currents rather than a non-zero value of the quark condensate. It is instructive 
to formally express Fq and the quark condensate by means of the euclidian functional integral 
of QCD in a finite box of size L. Integrating over the fermions first (giving them a mass m), 
the double limit lim m ^ h m L^oo of order parameters of SB^S becomes sensitive to the small 
eigenvalues X n (A) of the Dirac operator ip(A), averaged over all gluonic configurations. In 



particular, the crucial question is how dense does the spectrum become as L — > oo.p3|, |2|, [14 
If the average level spacing AA behaves as 1/L, then chiral symmetry remains unbroken. Both 
F and the quark condensate vanish. In order to have F ^ 0, it is necessary and sufficient that 
for L — > oo, AA behaves at least as 1/L 2 Fq [15|. It is well known that quarks start to condense 
only provided fL4| AA ~ 1/L 4 |(gg)|. The above discussion does not tell us what actually does 
happen in the QCD ground state. However, it demonstrates that the theoretical possibility of 
SB%S (i.e., Fq 7^ 0) without quark condensation might naturally arise within QCD. In fact, it is 
conceivable that small eigenvalues with AA ~ L~ 2 and AA ~ L~ 4 coexist and compete, giving 
rise to a marginal quark condensate, which then would be difficult to estimate in advance. 

At this point, one should recall that the alternative scenario of a small or vanishing 
condensate, though not contradicted by a single experimental fact, does not seem to be favoured 
by existing lattice simulations when they are extrapolated to the chiral limit. However, lattice 
regularisation is known to mistreat chiral symmetry, in one way or another, especially in the 
quenched approximation, in which the most significant results are obtained. (For a recent 
review, see Ref. |I6|.) Lattice results constitute another reason to attempt a direct experimental 
test of theoretical prejudices regarding the value of (qq). 

The quark condensate manifests itself exclusively through symmetry-breaking effects pro- 
portional to quark masses. The order parameter Bq (1.1) is one of the constants of the low- 
energy effective lagrangian, but it enters into physical observables multiplied by a quark mass. 
An example of a meaningful question to ask concerns, for instance, the magnitude of the 
renormalisation group invariant product (m u + rridjBo measured in units of the pion mass. 
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The smallness of most symmetry-breaking effects is the main reason why experimental tests 
of quark condensation are difficult and have not been performed so far. Actually, the only 
symmetry-breaking effects which are easily accessible experimentally are the Goldstone boson 
masses. Their expansion in powers of quark masses can be schematically written as (for more 
details, see fL7|, §) 

Ml + = {m u + m d )B + (m u + m d ) 2 A + ... (1.2) 
m k+ = ( m u + m s )B + (m u + m s ) 2 A + ... , 

where Fq A is a (suitably subtracted) massless QCD two-point function of scalar and pseu- 
doscalar quark densities. A is of the order of 1, typically Aq = 1 4- 5. (A similar expression 
holds for the r\ mass, except that at the quadratic level, a new constant appears, reflecting the 
axial U(l) anomaly.) While the overall convergence of the expansion (1.2) is controlled by the 
small parameter m quar k/A#, the relative importance of the first two terms is determined by the 
relative sizes of the quark mass and of the mass scale 

mo = ik (L3) 

set by the quark condensate. If m quar k <C m , the condensate term dominates the expansion 
of Goldstone boson masses. On the contrary, if m quar k ~ m , the first and second order terms 
in Eq. (1.2) are of comparable size. This is what distinguishes the standard large-condensate 
wisdom from the alternative, low-condensate, scenario. The standard case (mo ~ Ah) is the 
basis of the standard xPT, i.e. the expansion in powers of quark masses and external momenta 
in which one power of a quark mass is counted as two powers of the momenta: 

m quark = 0(p 2 ) . (1.4) 



The consequences of this scheme are well known. At leading order, pseudoscalar masses satisfy 
the Gell-Mann-Okubo (GMO) formula, and the ratio of strange to non-strange quark masses 
gets determined to be [TS[] 

r = ^ = 2 xf "1-26, m = -(m u + m d ). 1.5 
m M£ 2 

Historically, the emergence of the GMO formula at leading order was at the origin of the large B 
scenario 0. Today, the GMO formula can still be considered as an argument of plausibility, but 
not as an experimental test of the large-condensate hypothesis, for the following two reasons: i) 
The validity of the GMO relation does not imply the dominance of the condensate term in Eq. 
( |1.3| ) and it does not rule out the small condensate alternative, and ii) the standard is not 
yet able to predict the size of corrections to the leading order GMO formula; this would require 
an independent determination of the constant L§ |5fl from other symmetry-breaking effects. 
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On the other hand, if the quark condensate turned out to vanish, then the quark mass 
ratio (1.5) would become 



Mir 



6.3 . (1.6) 



Notice that even in this extreme case, the rj mass would remain unrelated to Mk and M n , 
because of the anomaly contribution. However, the small-condensate alternative is more general 
than this: it covers the whole range m <C Ah- For instance, if B ~ Fq ~ 90 MeV (i.e. 
— (qq) ~ Fq), then m could easily be as small as ~ 20 MeV. In this case, m quar k ~ m could 
happen even for u and d quarks, implying a value of the ratio (m u + itl^Bq/M^ significantly 
below one. In the generic case m <C Ah, the quark mass ratio r interpolates between the two 
extreme values (1.5) and (1.6). The standard xPT treats the terms m quar k/A# and m quar k/mo 
in a similar way: they are both considered small and of order 0(p 2 ). In the small-condensate 
alternative, the standard expansion of the (same) effective lagrangian has to be modified since, 
now, wiquark/^o is of order 1. This requires a modification of the chiral counting rule (1.4). Bq 
now becomes an expansion parameter that counts as a small quantity of order M n ||. The 
new counting rule reads 

m quark = O(p) , B = O(p) . (1.7) 

The corresponding expansion of C eS defines the so-called "generalized chiral perturbation the- 
ory" (G%PT) H H, [19], [l7i . Order by order, the G;\PT contains the standard expansion as 



a special case, because at each order it includes additional terms, which the standard xPT 
relegates to higher orders. Notice that in both schemes the symbol 0(p n ) has the same mean- 
ing: it represents a quantity with an order of magnitude (M n / A H ) n - The symbol 0(m" uark ), 
in contrast, has to be interpreted according to (1.4) or (1.7), depending on the scheme. More 
details on the structure of the G%PT may be found in Refs. || [19], [17|] . 

In order to test the strength of quark condensation, the constraints arising from the 
discussion of Goldstone boson masses have to be combined with experimental information on 
different symmetry-breaking effects. There are two particularly significant examples of such 
effects: 

i) Comparing the observed deviation from the Goldberger-Treiman relation in 7rN, KA and 
KYi channels one obtains a sum rule for the ratio r = m s /rh valid up to higher order 
corrections [^(J. The output of the sum rule is particularly sensitive to the precise value of 
the charged pion-nucleon coupling constant g n N, which is nowadays rather controversial. 



For the old value of Koch and Pietarinen ||21 1 , recently confirmed by the analysis of a 



new n — p charge-exchange scattering experiment , the value of r turns out to be 
significantly smaller than the standard value r ~ 26, by at least a factor of 2 [20]. For 
lower values of g n N the outcome is less conclusive. In any case, more accurate data as 
well as an estimate of higher order corrections are needed to transform this indication 
into a real measurement of m s /m. 
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ii) The second relevant symmetry-breaking effect concerns the 1=0 S-wave tttt scattering 
length. The standard x?T, based on the large-condensate hypothesis, predicts (at one- 
loop accuracy) Og = 0.20 ± 0.01 [23[ Q], while the current experimental value is 



a[] = 0.26 ± 0.05. One of the purposes of this work is to show that this one standard 
deviation effect finds a natural explanation within the low-condensate alternative. It 
has been known for some time || that already at tree level Weinberg's prediction [EE 



a[] = 0.16 (which implicitly uses the large-condensate hypothesis) can gradually increase 
towards a value a(j = 0.27 provided the ratio {m u + md)Bo/M% is allowed to decrease 
from 1 to 0. Here, this statement will be made more quantitative including the one-loop 



26 1 and the two-loop G^PT corrections. In view of forthcoming precise low energy ixix 
scattering data |TI| [Ll| , our analysis might imply a realistic possibility to disentangle the 
large- and low-condensate alternatives experimentally. 

The paper is organized as follows. We first describe the relevant features of the effective 
action at order 0(p 4 ) in the generalized case and obtain the corresponding expression for the 
ixix amplitude in terms of four independent combinations of the low-energy constants (Section 
2). A preliminary version of this calculation has already been published [p6 ]. In Section 3, we 



show how this result can be rederived and easily extended to two-loop accuracy using a slightly 
different approach. At two-loop order, the tttt amplitude involves two additional constants 
which are not fixed by chiral symmetry. We determine four out of this total of six constants 
from sum rules and tttt production data at higher energies in Section 4, where the dependence 
of low-energy phase shifts and threshold parameters on the remaining two parameters is also 
discussed. The possibilities of exploiting future precise data on low-energy tttt scattering in 
order to distinguish between the small and large condensate alternatives are investigated in 
Section 5. A summary and concluding remarks are presented in Section 6, and miscellaneous 
formulae and auxiliary results have been collected in four Appendices. 

Finally, it should be stressed that although the paper is motivated mainly by questions 
raised by the chiral structure of the QCD vacuum, this work contains new material concerning 
the phenomenology of low-energy ixix scattering which might be of interest for its own sake. 
The reader mostly interested in these aspects can directly read Sections 3 and 4 as well as 
appendices B, C and D, which form a self-contained core. 

2 The tttt Amplitude to Order 0(p 4 ) in G%PT 



In this section, we briefly discuss the structure of the generating functional for the effective 
theory in the generalized framework up to order 0(p 4 ). Then we derive the low-energy expan- 
sion, to this order, of the amplitude A(s|t, u) which describes ixix scattering in the absence of 
isospin breaking (i.e. we shall take m u = rrid = m and we shall neglect electromagnetic correc- 
tions). We work within S'?7(3)l x SU(3)r, although ixix scattering is merely an SU(2) L x SU(2) R 
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problem. The reason for this choice should be obvious from the preceding discussion on Gold- 
stone boson masses. Working within the S77(3)l x £77(3)r chiral expansion offers the opportu- 
nity to make use of relationships between constants which characterise low-energy tttt scattering 
and QCD parameters (such as m s /m) or observables (such as Fx/F n , Mx/M n ). 

The strength of the quark-antiquark condensate in the QCD vacuum does not affect the 
structure of the effective lagrangian £ e£f . The latter is merely determined by chiral symmetry 



27 1 and the transformation properties of the chiral symmetry breaking mass term of QCD. 



Thus, £ consists of an infinite tower of invariants 



•.eff 



(M) 



where £<k,i) contains k powers of covariant derivatives of the Goldstone boson fields, and I 
powers of scalar or pseudoscalar sources. In the chiral limit, £(k,i) vanishes like the k-th power 
of external momenta p and the /-th power of the quark masses, 

~ PHuark • (2-2) 

For sufficiently small quark masses, such that both m quar k <C A H and m quar k <C m Q = B Q /2A , 
one has £ = 0(p k+21 ). In this case, the double expansion (|2.1| ) can be reorganized as 

£ eff = £ (2) + £ (4) + £ (6) + . . . ^ (2 _ 3) 

where [g, § 

£ [d) = E Am- (2-4) 

k+2l=d 

This expansion defines the standard x?T. On the other hand, if for the actual values of the 
quark masses one has m quark ~ m Q A H , then £(h,i) = 0(p k+l ), and this new counting yields 
a different expansion of the same effective lagrangian ( |2.1| ) : 



£ cS = £^ + £(3) + £W + £(5) + £(6) + . . . ) (2 . 5 ) 

with now || 

£«= E B ^(k,i)- (2-6) 

k+l+n=d 

Although Eqs. ( p.3| ) and ( |2.5| ) sum up the same effective lagrangian £ cS to all orders, their 
truncations at any finite order may differ. 

The structure of the effective action Z in the generalized case has been discussed in 
up to order 0(p 4 ) in the framework of the S'?7(3)l x SU(3)r chiral expansion. It is given by 

Z cff = / d 4 x [£^ + £^ + m} + z[f oop . (2.7) 
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Whereas C^ n \ n = 2,3,4, gives the tree level contributions at order 0(p n ), z[f oop contains the 
contributions from one-loop graphs with an arbitrary number of vertices from only. Thus, 
the leading order of the generalized expansion is described by 2 \ which was first given in ||: 

£~( 2 ) = \f 2 {{D,U + D»U) +2B (U + x + X + U) + 

+A ((U + xY + {X + Uf) + Z S (U + X + X + U) 2 (2.8) 
+ Z?(U + x-X + U) 2 + H (x + x)} ■ 

The notation is as in Refs. [Q, except for the consistent removal of the factor 2B from x, 
the parameter that collects the scalar and pseudoscalar sources, 

X — s + ip = M H , M = diag(m u , m d , m s ) . (2.9) 

In G^PT, the next-to-leading-order corrections are of order 0(p 3 ), and still occur at tree level 
only. They are embodied in CS 2 "' = £(2,1) + £(0,3), which reads || [T7| 

£ (3) = \Fq {t(D,U + D»U(x + U + U + X )) +i(D,U + D»U){ X + U + U+ X ) 
+Pi((x + U) 3 + (U + X f) + P2((X + U + u + x )x + x) 
+ P3 (x + u - u+x)((x + u) 2 - {U+x?) (2.10) 

+P4((X + U) 2 + (U+x?) (x + u + U + x) 

+ps{x + x)(x + u + u + x) 

+ p 6 ( X + U - u + x ) 2 (x + u + u + x ) + p 7 (x + u + U + x) 3 } ■ 

The tree-level contributions at order 0(p 4 ) are contained in 1 

£ (4) = £( 4j0 ) + £( 2 ,2) + £(0,4) + -Bo£(0,2) + 5 o£(2,l) + A)£(0,3) ■ ( 2 -H) 

The part without explicit chiral symmetry breaking, £(4,0), is described by the same low-energy 
constants L\, L 2 , L3, L 9 and L 10 as in the standard case |5[]. The part with two powers of 
momenta and two powers of quark masses is given by: 

£(2,2) = \f 2 ( [ a 1 (d^u + d^u(x + x + u + xx + u)) 

+A 2 (D ll U + U X + D»UU + x) 

+A 3 (D,U + U(x + D»X ~ D»x + X) + D,UU + {xD^x + ~ D»XX + )) 

+A 4 (D,U + D» U) ( X + x) 

x At order 0(p 4 ), terms with odd intrinsic parity coming from the Wess-Zumino term are also present in the 
effective lagrangian, but since they start to contribute to A(s\t, u) only at order 0(p 8 ), we do not mention them 
further here. 
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+B 1 (D^U + D»U(x + U X + U + U + X U + x)) 
+B 2 {D,U + X D»U + x + X + D,U X + D»U) 
+B i {D,U + D»U)( X + Ux + U + U + X U + x) 

+C?(D fl U X + + xD»U+)(D»U X + + XD»U + ) (2.12) 
+D s (D,U + D»U(x + U + + t/ + X> + •••}. 

Finally, the tree- level contributions which behave as 0(m 4 uark ) in the chiral limit are contained 
in £(o,4), which reads 

£(0,4) = ^o 2 {^i((x + ^) 4 + (^ + x) 4 ) 

+e 2 (x + x(x + u x + u + u + x u + x)) 

+E 3 ( X + xU + XX + U) 

+Fi s ( X + Ux + U + U + X U + x) 2 
+F 2 s ((x + U) 3 + (U + X ?)(X + U + U+x) 

+Fi( x + x(x + u + U+x)) (x + u + U + X ) 

+F 4 s ((x + U) 2 + (U + X ) 2 }(x + X} 

+Fi s ((x + U) 2 + (U + X ) 2 )(X + U + U + X ) 2 (2.13) 

+F 6 ss (x + x)(x + U + U+x) 2 + ---} ■ 

For £(2,2) and £(o,4), we have shown explicitly only those terms which will actually contribute 
to A(s\t, u). For the complete list of counterterms which enter £(2,2) and £(0,4), we refer the 
reader to ||28|| . Notice also that in the standard framework all these contributions would count 
as order 0(p 6 ) and order 0(p 8 ), respectively Q. 

As far as the mr amplitude is concerned, the 0(p 4 ) loop corrections involve only graphs 
with a single or two vertices from C^: 

£(4) _ 2 (4) 4- # (4) -I (2^4-) 

^lloop ~~ ^tadpole * ^unitarity ' K*" 1 *) 

The divergent parts of these one loop graphs have been subtracted at a scale \i in the same 
dimensional renormalisation scheme as described in 0, M. Accordingly, the low energy constants 
of £(4,0), £(2,2), and £(0,4) stand for the renormalised quantities, with an explicit logarithmic 
scale dependence (X(fi) denotes generically any of these renormalised low-energy constants) 

XQi) = X(fx') + • ln(////i) . (2.15) 

At order 0(p 4 ), the low-energy constants of and £ ( - 3 ^ also need to be renormalised. The 
corresponding counterterms, however, are of order O(Bq) and O(B ), respectively, and they 
2 The standard 0(p 6 ) effective lagrangian = £(6,0) + £(4.i) + £(2,2) + £(0,3) has been worked out in Ref. 
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are gathered in the three last terms of Eq. ( |2.11| ): in G%PT, renormalisation proceeds order by 
order in the expansion in powers of Bq. Alternatively, one may think of Eqs. (|2.8|) and fl2.10|) 
as standing for the combinations + BqC'^ 2 - ) and +B £', 21 s +B C'^ 3 ^, respectively, with 
the corresponding low-energy constants representing the renormalised quantities. The full list 
of /3-function coefficients Fx can be found in Ref. [p8fl , where the explicit expressions of tadpole 
and of ^unitarity are a l so displayed. 

Computing the 0(p 4 ) expression of ^4(s|t, u) is then a straightforward exercise in field 
theory, and we shall merely quote the result, expressed in terms of the loop functions J t P q and 
MpQ defined as in g §: 

B r / 4 \ /If 2 



+ ±(2L[ + L 3 )(s-2M 2 ) 2 + irU 



F 4 



(t - 2M 2 ) 2 + (u - 2M 2 ) 2 
+ Wa JUs)[A[s - -f + 10m 2 (A + 2Z*)} 2 -[s- 2M 2 - 8m 2 (A + 2Z S )] 



+ 
+ 
+ 
+ 



4F 4 
1 

8F* 
1 



Jkk{ s 



t - 2Ml - 8m 2 (A + 2Z$ 



+ 



1 



4Fi 



M-2A4 2 -8m 2 (A + 2Zo' 



ls/ .4 J U S 



2Fi 



s + 8m 2 (l + r)(A + 2Z s 

M 2 + 4m 2 (1 + 2r)(A + 2Z%) + 8m 2 (1 - r)(A + 2Z? 
s - u)t[2M^{t) + M r KK (t)\ + (s- t)u[2M^(u) + M r KK (u)]\ . 



(2.16) 



The coefficients a r and j3 r denote two combinations of the (renormalised, scale dependent) 
low-energy constants of + + 



(3 r 



1 + 2m£ + 4m£ 



+2m" 



3A 2 + 2A 3 + AB X + 2B 2 + 8B 4 + 4Cf + 2(2 + r)D 5 
-2e 2 -4(2 + r)e 2 -2(2 + r)(4 + r)^1 , 



(2.17) 



and 



p2 

— M 2 a r 
-^0 



2mS + 16m 2 A + 4(8 + r)m 2 Z$ - 4mM^ + 2f) 



+ m d 



81pi +p 2 + 2(82 + 16r + r 2 )p 4 + (2 + r 2 )p 5 + 12(2 + r)(14 + r)p 7 



- 8m 2 M 2 \2B 1 - 2B 2 + 4B 4 - A 3 - 4Cf + (4 + r)L> 5 

+ 2m 4 [A8A 3 (A + 2Z^) + 128E 1 + 8E 2 

+32(8 + r 2 )if + (272 + 81r + r 3 )F 2 5 + (16 + r + r 3 )F 3 s + 8(2 + r 2 )Ff 
+4(144 + 82r + 16r 2 + r 3 )F 5 5S + 2(16 + 2r + 8r 2 + r 3 )F 6 5Sl 
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^[3M 2 + 116m 2 (A + 2Z S )] - 2^[M 2 + 4(7 + r)m 2 (A + 2Z 5 )] (2.18) 
- ^ V [M 2 + 4(7 + 2r)m 2 (A + 2Z 5 ) + 8(4 - r)m 2 (A + 2Z?)] , 

with iip = Mp/ (32n 2 F 2 ) \n(Mp / fi 2 ), P = 7i,K,r]. The deviation of the pion decay constant 
F n from its counterpart in the chiral limit F is given, at order 0(p 4 ), in Appendix A, together 
with the corresponding expansion of the pion mass. 

As mentioned above, both a r and (3 r depend of the same subtraction scale fi in terms of 
which the renormalised loop functions J PQ and M PQ were defined. Using the results of 
Ref. H§, we find [] 



M 2 r 
aTO*) = a V) + ^l«'{[l + (l + ^[l+3(l + ^ 



(4ttF 7 

+hl + (1 + 2r)e + 2^^] 2 + [1 + 22e + 33e 2 ] ) , (2.19) 



and 



M 2 

= /7(//) + — Ih(ai7ai) [5 + (11 + r)e] , (2.20) 



where we have introduced 



and 



— (A + 2Z 5 ), (2.21) 



A G mo = A{ ™ M T S? (Ao + 2Z P ) . (2.22) 

7T 

From these formulae, one verifies that A(s\t, u) is indeed independent of the subtraction scale 
/i, thus providing us at the same time with a nontrivial check of our calculation. 

Using the 0(p 2 ) expressions for the pseudoscalar masses derived from |20], ||, |9|, |28[| , 
one obtains 

? = 2 3^t( 1 + 2 t 1 + 0(m quark )] , 

(2.23) 

A GMO = (3M 2 -4M| + M 2 )/M 2 [1 + 0(m quark )] , 

with r 2 = 2M\jM 2 — 1 ~ 26, £ = Zq /A , which is expected to be small due to the Zweig rule, 
and r = m s /m. In G^PT, the latter quark mass ratio remains a free parameter. 

On the other hand, in the standard case, r is determined order by order in the chiral 
expansion. For instance, at leading order, r = r 2 (and A GMO = 0), while at order 0(p 4 ), one 



3 Up to order 0(p 3 ), a r and (3 r are scale independent and coincide with the parameters a^^ and (3 V7! which 
enter the 0(p 5 ) expression of the 77 — > 7r°7r° amplitude in generalised x?T p9|. 
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obtains |5| 



' St 



'Ml 



A 



M 



(2.24) 



A 



v = — + ^ + — 2 {M 2 K - M 2 )(2L^(/i) - I7M) 



Furthermore, the contributions from £(0,3), £(2,2) and £(0,4) are relegated beyond order 0(p 4 ). 
The remaining constants A , Zq , Zq, £ and £ are related to the low-energy constants Lj of 
introduced by Gasser and Leutwyler 0, by (both sides of these equations refer to the 
renormalised quantities, defined at the same scale /i): 



w 2 



8B 



F 2 



16B 2 ' 



T T 



8B 



F 7I' 



(2.25) 



IQB 2 



^0 > ^8 



Thus, upon replacing r by its standard expression r st above, and upon keeping only those 
contributions that are of order 0(p 4 ) in the standard counting, one recovers from Eqs. (|2.16|) , 
( |A.1| ), ( |A.4| ) and ( |2.25| ), as a special case, the expression of the one-loop nit amplitude in the 
SU(3)l x S77(3)r framework of the standard chiral expansion |30fl : 



A(s\t, u)\ 



St 



F 2 



M 2 



+ + L 3 )(s - 2M 2 ) 2 + ±1% [(t - 2M 2 ) 2 + (u- 2M 2 ) 2 



+ 
+ 
+ 
+ 



2F 4 
1 

4^ 
1 

1 

2Fl 



s 2 -Mi 



(t-2M 2 YJUt) + (u-2M 2 yjUu) 

1 



s 2 JkM + 



18Fi 



(2.26) 



{{s - u)t [2ACW + M r KK {t)\ +(s- t)u [2MIM + M r KK (u)}} 



with 



St 



1 + ^(2^4+^5) , 



(2.27) 



a 



St 



16M 2 48M 2 



(2.28) 
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The scale dependences of a r st and /3* t are readily obtained, either from Eqs. ( |2.19| ) and ( |2.20| ), 
upon taking e and Aqmo equal to zero, or directly from the scale dependences of the L^s as 
given in Ref. 0. Both ways lead to the same result: 

7M 2 

(2.29) 

= /W) + (^fp 111 ^//^) • 

Starting from Eq. ( |2.26| ) above, it is then straightforward to further reduce this expression 
to the SU(2) L x SU(2) R case, and thus to recover the result, now expressed in terms of the 
constants k, i = 1,2,3,4, first obtained by Gasser and Leutwyler [^, f|. We shall not pursue 
this matter for the moment. 

3 The Use of Analyticity, Unitarity and 
Crossing Symmetry in xPT 



In the previous section, we have obtained the tttt amplitude to order 0(p 4 ) in the gener- 
alised chiral expansion, and have compared it to its standard counterpart. At this stage, two 
features are worth being kept in mind: First, in both cases, the amplitude depends on only four 
independent combinations of low-energy constants of the respective 0(p 4 ) truncations of the 
(same) effective lagrangian. Secondly, the essential difference between the two cases arises from 
the fact that the quark mass ratio r = m s /m remains a free parameter at each order of the 
generalised expansion. Upon fixing this ratio to its value predicted by the standard expansion, 
one recovers the nir amplitude in the latter framework as a special case. These observations 
lead to the expectation that the form of the titt amplitude might also be obtained beyond the 
0(p 4 ) approximation independently of any prejudice concerning the size of the quark-antiquark 
condensate. In this section we shall establish the expression of the 7m amplitude up to correc- 
tions of order 0(p 8 ). More phenomenological aspects, such as how the actual value of m s /m 
might actually be inferred from low-energy nn data will be the subject of subsequent sections. 

It has indeed been shown || that the Goldstone nature of the pion, combined with 
analyticity, crossing symmetry and unitarity of the S-matrix, determines the low-energy Tin 
scattering amplitude up to six arbitrary constants (not fixed by these general requirements) 
and up to corrections of the order 0([p/A#] 8 ). (As usual, Ah denotes the mass scale at 
which particles other than Goldstone bosons are created, Ah ~ lGeV .) The G%PT one- 
loop 0(p 4 ) amplitude obtained in the preceding section has to be a particular case of this 
general six-parametrical formula. Moreover, since the general low energy representation of 
the amplitude extends up to and including order 0(p 6 ) (actually 0(p 7 ) in the generalised 
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counting), it must in fact contain the whole two-loop amplitude both in the standard and 
in the generalised settings. In this section, we shall first rederive the 0(p 4 ) amplitude using 
dispersion relations, and subsequently these same dispersive methods will be used to work out 
the two-loop amplitude explicitly. 

The starting point is the following Theorem, proved in Ref. |§: For pion momenta 
p <C Ah, the 7T7T scattering amplitude can be expressed as 

^A(s\t,u) = I [W (s) - W 2 (s)} 

+ ^[3(s-u)W 1 (t) + W 2 (t)\ (3.1) 
+ ^[3(s~t)W 1 (u) + W 2 (u)] + 0([p/A H } s ) . 

The three functions W a (s), a = 0,1,2, are analytic except for a cut singularity at s > 4M%. 
Their discontinuities across this cut are given by 

lmW a (s) = {s-AMl)- e «\mf a {s)6(s-AMl) , (3.2) 

where 

£o = e 2 = , ei = 1 , (3.3) 

and with f a (s), a = 0, 1,2, denoting the three lowest partial-wave amplitudes //(s) (the 1 = 
and 1 = 2 S- waves, and the P-wave), 

fo(s) = / °(s) , fi(s) = /i(s) , f 2 (s) = / 2 (^) • (3-4) 

(Notation, normalisation and other conventions are summarised in Appendix A of Ref. 0.) 

The proof of this Theorem has been given in Ref. |5|] with all details and it will not be 
reproduced here. Let us just recall that the main ingredient of this proof is the suppression of 
inelasticities and of higher partial waves in the chiral limit p — > , M n — > 0, p/M n fixed: Below 
the threshold of production of non-Goldstone particles, one can write 



Im /a ( s ) = ^^_ ^| /a(s) |2 + o([p/A„] 8 ) . (3.5) 

Similarly, the absorptive parts of higher partial waves are suppressed as 

lmf!(s) = 0([p/A H } 8 ) , l>2. (3.6) 

As discussed in Ref. 0, the properties ( |3.5|) and ( |3.6| ) are rigourous consequences of the fact 
that the whole amplitude A(s|t, u) behaves dominantly as 0(p 2 ) in the chiral limit, combined 
with analyticity, unitarity, crossing symmetry and chiral counting. At the same time, these two 
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equations explain why it is difficult to extend the representation (|3.1|) beyond two loops: At 
order 0(p 8 ), inelasticities and higher partial waves set in. 

Given Im/ a (s), then the four functions W a (s), and consequently the whole amplitude 
A(s|t, u), are determined up to a polynomial. The maximum degree of this polynomial may 
be fixed according to the chiral order 0(p 8 ) of the neglected contributions: Without loss of 
generality, the functions W a (s) can thus be further specified by the asymptotic conditions 

lim s £a ~ A W a (s) = . (3.7) 

|s|^oo 



It then follows that for a given set of imaginary parts Im f a (s) consistent with the bounds (3.7), 
the amplitude A(s\t, u) is determined up to six arbitrary constants. 

It is worth stressing that in spite of Eq. ( |3.2|) , the functions W a (s) do not coincide with the 
corresponding partial wave amplitudes f a {s). The latter are obtained by the standard partial 
wave projections of Eq. (|3.1|) , 



= ^C,I^ fl ( 1 + 7^)™ t ') + - 4 < t ' s - 1 '» + ^ i »> 

/ ' (s) = m L ; -, db p < i 1 + j^m) (A(i|s < »> - A ^ i)} (3 - 8) 

where the integrals are taken with s kept fixed and with u reexpressed as u = 4M% — s — t. The 
functions W a (s) have only a right-hand cut for s > 4M^, whereas the partial wave amplitudes 
f a (s) Q3.8Q do have both a right-hand cut (s > &M%) and a left-hand cut (s < 0), with the 
correct discontinuities, as dictated by unitarity and crossing symmetry. Actually, the partial 
waves with I > 2 obtained from ( ^.1[ ) are real, in agreement with Eq. ( |3.6| ). The formula ( |3.1| ) 
may be viewed as the most general solution of the constraints imposed by analyticity, unitarity 
and crossing symmetry up to and including the chiral order 0(p 6 ) (0(p 7 ) in the generalised 
case). In practice, the representation ( ^.1| ) should be a good approximation within a range of 
energies where the inelasticities in ( |3.5| ) and the contributions of higher partial waves in ( |3.6| ) 
remain small. Let us recall that the inelasticities, though in principle present above 540 MeV, 
are actually observed to be very small even up to 1 GeV. The same remark applies to the higher 
partial waves, I > 2. 



3.1 First Iteration of the Unitarity Condition 

The dispersion relation approach to is based on the iteration of the unitarity con- 
dition. Here, the latter will be considered in the elastic form fl3.5p , i.e. below the KK and rjrj 
thresholds. This does not mean that we have to renounce the advantages of an SU (3)l x S77(3)r 
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analysis of tttt scattering, but rather that the K and r\ loops will be expanded in powers of 
s/{2Mk) 2 and s/(2M r) ) 2 , respectively. Since in practice 2Mk and 2M V are of the order of A#, 
this should not affect the accuracy of our results at low energies E <C A H . 

Chiral symmetry implies that the low-energy expansion of the irn amplitude starts at 
order 0(p 2 ). Analyticity and crossing symmetry then restrict its form to 

A(s\t, U ) = a^ + P 2 (s- ~M 2 ) + 0(p 4 ) , (3.9) 

with a and f3 two constants whose dominant behaviour in the chiral limit is of order 0(1). 
In G^PT, where odd chiral orders also occur, the constants a and (3 in Eq. ( pT9| ) consist 
of both 0(1) and 0(p) contributions, the latter being proportional to the first power of the 
quark masses. The Goldstone nature of the pions together with the general properties of their 
scattering amplitude, by themselves, do not fix the constants a and (3. On the other hand, the 
values of these constants are actually at the center of our interest. They reflect the mechanism 
of spontaneous chiral symmetry breakdown in QCD and the corresponding chiral structure of 
the ground state: For instance, in the limit of a large quark condensate (B ~ A#, the special 
case of standard x?T), one has a — 1, /3 — 1, whereas in the opposite limit of a vanishing 
condensate, a = 4 and f3 = 1 at order 0(p 2 ). It is important that the dispersive approach to 
%PT and, in particular, the iteration of the unitarity condition, only assumes that a and f3 are 
of the order 0(1) in the chiral counting, but otherwise proceeds independently of any particular 
values of these constants. 

Equation ( |3.9|) implies for the three lowest partial wave f a (s) , a = 0, 1, 2, 

Re f a (s) = Ms) + 0(p 4 ) , Im/ a (s) = 0(p 4 ) , (3.10) 



where 



L_ | 6 /?( S -^M 2 ) + 5aM 2 } 



96ttF 2 
1 



Ms) = 9^/^- 4M *) (3- 11 ) 



Ms) = 9 ^{-^( s -^ 2 ) + 2 ^ 2 } • 

Inserting this information into the unitarity condition (|3.5|), one obtains 



lmf a (s) = ^-—^l^is)] 2 + 0(p 6 ) , (3.12) 

which are the discontinuities of the functions W a (s) up to and including the order 0(p 4 ~ 2e ") 
(0(p 5 ~ 2ea ) in G^PT). Due to the polynomial character of <f a {s), Eqs. (|3.2j ) and (|3.12| ) allow 
one to express the functions W a (s) as 

W a (s) = 16vr (s - AMlY £a [M s )} 2 J(s) + polynomial + 0(p 6 - 2£a ) , (3.13) 
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where J(s) denotes the dispersion integral 



J(s) 



16tt 2 



dx 



I x 



AMI 



AMI X X — S 



X 



(3.14) 



The reader will easily recognise in the analytic function J(s) the standard scalar-loop integral 
J^ n (s) of the previous Section, subtracted at zero momentum transfer, J^{s) = J^ n (s) — 7^(0), 
(see Appendix C). The polynomial in Eq. (|3.13|) is at most of degree 2 — e a in s. The coefficients 
of higher powers of s should be well behaved in the limit M n — > , for s fixed, and consequently, 
such higher powers can be absorbed into the neglected 0(p 6_2ea ) remainders. Notice that in 
the case of Eq. (|3.13| ), the asymptotic bounds ( |3.7|) are not only satisfied, but they are not even 
saturated (J(s) grows as Ins ). Returning to the formula ( p.l|) , the polynomial contributions 
in W a (s) beyond the leading contributions ( |3.9| ) can be collected into 

4 
3 
A 2 



5A(s\t, u) 



-m: 



2\2 



(t - 2M^Y + (u- 2M 



2\2 



K> P4 

Here, in the chiral limit, the dominant behaviour of 8a, 8/3, Ai and A 2 is 



(3.15) 



5a = 0(p 2 ) , 5(3 = 0(p 2 ) , A li2 = 0(1) . 



(3.16) 



It is convenient to absorb 5a and 5(3 into the constants a and (3 which characterise the lead- 
ing order contribution (|3.9| ). The constants a and (3 so redefined by including higher-order 
contributions will keep their original names. Notice that in the expression ( pj,13| ) for W a (s), 
the redefinitions a — > a — 5a, (3 —>■ (3 — 5(3 induce terms of higher chiral orders, which can 
be reabsorbed into the neglected 0(p 6_2£a ) remainders. Hence, the whole result of the first 
iteration of the unitarity condition can be expressed in terms of four parameters, a, (3, \\ and 
A 2 : 

Ml 



A(s\t,u) = 



+ 
+ 



+ a 



s - 2M 2 ) 2 + 



A2 

F A 



(t - 2M 2 ) 2 + (u- 2M 



2\2 



J(a,8)(s\t, U)+0( P 6 /A G H ) 



(3.17) 



where 



Ji 



s\t, u) 




/3( S -^M 2 ) + ^«M 2 ] 2 - 
W-\Ml)-\aMl 2 



- -Ml) 



zOiMt 



P{s-\Ml)-\aM^\j{s) 



+ (3 2 {s-u){t-AMl))J(t) (3.18) 



+ (3 2 {s-t){u-AMl)\j{u) 
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3.2 Comparing Perturbative 0(p 4 ) and Dispersive Formulae 



The parameters a, (3, Ai and A2 in Eqs. (|3.17 ) and (|3.18| ) have their own expansions in 



powers of quark masses and they involve chiral logarithms. One has 

«=X> (n) , (3=jzP {n \ A 1>a = (3.19) 

n=0 n=0 n=0 

where a^ 1 ' denotes the 0(p n ) component of a and likewise for the other constants. In the 
standard x?T, the odd orders are absent, = (3^ = 1, and the formula ( |3.17|) contains 0(p 2 ) 



and 0(p 4 ) parts only. In generalised xPT, the odd orders show up, (3^ = 1, and 1 < cr°) < 4, 
whereas the amplitude ( |3.17| ) involves 0(p 2 ), 0(p 3 ), 0(p 4 ) and 0(p 5 ) contributions before it 
reaches the (neglected) order 0(p 6 ). Detailed discussion of the values of the constants a, j3, 
A12 will be given in Sections 4 and 5 . 

It is now easy to show that the 0(p 4 ) result ( |2.16|) obtained by a direct G^PT evaluation 
of one- loop graphs, starting from the Lagrangian + C^' + is indeed of the form of Eqs. 
(|3.17|) and (|3.18|) . To see this, one should first reexpress the scale dependent renormalised pion 
loop integrals Jl^s) and M^(s) in terms of the scale independent function J(s): 




MlAs) = ^ (s - 4M 2 ) J(s) 



(3.20) 




Next, the K and r\ loop contributions J T KK {s\ M r KK (s) and J^ v (s) have to be expanded in 
powers of sj (2Mk) 2 and s/ (2M V ) 2 , respectively. Keeping only terms that contribute up to and 
including the order 0(p 4 ), one has 



1 / M 2 \ 
J r p P (s) = f In + 1 J + 0( S /4M|) , P = K, V 



1 / M 2 \ 
M kk(s) = (in -f + lj + 0( S /4M£) 



(3.21) 



Let us stress again at this point that it is perfectly possible to recover the full K and r\ loop 
contributions within the dispersive approach: for this, it is sufficient to include the KK and rjrj 
intermediate states into the unitarity condition. Here we do not present this extended analysis 
for reasons of simplicity, and because performing the expansions ( |3.21| ) in Eq. ( |2.16| ) should be 
a fairly good approximation at low energies. 

Finally, in order to compare the two expressions, all contributions beyond the order 0(p 4 ) 
in Eqs. ( |3.17| ) and ( |3.18[ ), i.e. all 0(p 5 ) contributions, should be dropped. In particular, only the 
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0(p 4 ) part of the loop function J( a ^)(s\t, u) should be maintained. This amounts to replacing 

J(a,f3)(s\t, u) by J(«(o) jj8 (o))(s|t, u), where 



a 



(o) 



l + 3e, (3 



(0) 



(3.22) 



are the leading parts of the constants a and /?, as can be inferred from Eqs. Q2.17|) , (|2.18 
Then one finds that Eqs. (|2.16|) and (|3.17| ) are indeed the same, upon making the following 
identifications (e and A GMO were defined in Eqs. ( |2.21| ) and ( |2.22j )) 



a: 



c/(/i) 



Ml 



32tt 2 F 2 
Ml ' 



96tt 2 F 2 
Ml 



Ml 
In -3- +1 



+ 1J [1 + (1 + r)e] [1 + 3(1 + r)e] 

-i 2 



327T 2 F 2 i 111 



Ag 

2 



= - 



Ml 



l6n 2 Fl 



In 



/i 2 



Ml 



32vr 2 F 2 



In 



Mi 



A' 



l + (l + 2r)?+2 
l + 22e + 33e 2 " 

+ 1 j [2 + 5e] 
1) [l + (l + r)e] , 



A 



GMO 



1 



(3.23) 



(3.24) 



and 



Ai = Ai' 



(0) 



A 2 



.(0) 



4(2LU/i) + L 3 ) 



48vr 2 



48vr 2 



\ M 2 1 Mf 
In — f + - In ,v 



1 M| 35 

- In — £ H 

8 /x 2 24 



23 



(3.25) 



/i 



J + 24 



Upon using Eqs. ( |2.19|) and ( |2.20| ) and the scale dependences of the renormalised constants L\ 
and L r 2 , it is straightforward to ascertain the scale independence of a, (3, \\ and A2. 

At the end of Section 2 it was shown how the standard SU (3)l x SU (3)r 0(p 4 ) amplitude 
appears as a particular case of Eq. ( |2.16| ): It may be expressed in terms of the standard 0(p 4 ) 
constants Li,... L 8 . The expressions ( |3.25|) for Ai and for A2 remain unchanged, while one finds 
that a and (3 become (a^ = /3, 



i 0) 



M 2 M 2 
16^(2Ll + L r 8 ) + 4&^(2Ll + Ll 



1 Ml 



Ml 



< hi + In 



M 



a' 



tit 



32vr 2 F 2 
M 2 

8-^(2^ + L 



1 , mi r 

+ - In —3- + - 

3 /i 2 3 



(3.26) 



32vr 2 F 2 



M 2 f , M 2 M 2 



(3.27) 
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If we use for the constants the values currently given in the literature |!T| and take = 
92.4 MeV, M n = 139.57 MeV, we obtain, for the corresponding scale independent constants a, 
(3, Ai and A 2 the following values at order 0(p 4 ) in the standard case: 



At 

Ai,st 

^2,st 



1.04 ±0.15 
1.08 ±0.03 

(-6.4 ±6.8) x 10" 
(10.8 ± 1.2) x 10" 



(3.28) 



Similarly, the standard SU(2) L x S77(2) R 0(p 4 ) amplitude originally given by Gasser and 
Leutwyler [^3|, f| in terms of the four scale independent parameters h, I2, h and U can be put 
into the form ( |3.17|) and (|3.18| ), with the identifications 

«gl = 1 + 3^ |f (4/4 - 3/3 - l) , (3.29) 



0< 



GL 



and 



Ai 



,GL 



48vr 2 



1 3 



A 



2.GL 



48vr 2 



Using the numerical values recently updated by Gasser in f32 | 

h = -2.15 ±4.29 , l 2 = 5.84 ± 1.72 , l 3 = 2.9 ± 2.4 , l 4 = 4.55 ± .29 



(3.30) 



(3.31) 



(3.32) 



one obtains 



Al,GL 
A2,GL 



1.06 ±0.06 
1.103 ±0.008 

(-7.35 ±9.06) x 10~ 3 
(10.57 ±3.63) x 10^ 3 . 



(3.33) 



It thus appears clearly that the situation of a large quark condensate, i.e. the standard version 
of xPT, is characterised by the values of the parameters a and j3 remaining close to unity. 



3.3 The Two-Loop tttt Amplitude 

We now proceed to the second iteration of the unitarity condition ( |3.5|) . The starting 
point is the result of the previous step, i.e. the one-loop amplitude ( [3.171) . First, we use Eqs. 
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( |3.8| ) to calculate the S- and P-wave projections of the amplitude ( |3.17| ). All the integrals are 
elementary and the result can be represented as 

Re f a (s) = <p a (s) + ip a (s) + 0(p 6 ) , (3.34) 

where ip a {s) are the polynomials and the new contributions ij) a (s), which will be given 

shortly, are dominantly of order 0(p 4 ). The unitarity condition ( |3.5|) now allows us to push 
the knowledge of Im/ a (s) one step further: 



lmfa(s) = f A s M " {\Ms)\ 2 + 2^ a (s)Ms)} + Q(P 8 ) • (3.35) 

It is convenient to express the functions ip a {s) arising from the partial wave projections in the 
following way (s > 4M 2 ): 

/If 4 / s 4 

Ms) = E tf°00 Us) ■ (3.36) 

Here, £i ( s ) are polynomials in s of at most second degree, with coefficients given in terms of 
a, j3, Ai and A2. They are listed in Appendix B. k n (s) represent the following set of elementary 
functions 



, , x 1 Is -AMI , / x 1 / x 

107T V S 07T 

1 / 4/W 2 \ 3 M 2 

fe( a ) = - 1 * , **(*) = — ' L 2 (s) , (3.37) 

87T V S y 167T J s ( s _ 4M 2) 



, / x 1 AC s / x Mi 9 , s 

Ms) = is? 1 1 + i(s) + ^ 



where 



L(s) = In \ , s > 4M 2 . (3.38) 

i + v /rnM' 

Notice that in the chiral limit (s —>■ 0, M 2 — > 0, s/M% fixed) all functions fc n (s) are of order 

0(1). Furthermore, they are all asymptotically bounded by Ins. Eqs. ( |3.34|) and ( |3.35| ) now 
lead to the following expressions for lmW a (s) 

lmW a (s) = j2 w<?\s) k n (s) + 0(p 8 - 2£ «) , (3.39) 

n=0 



where 



{9 /if 4 1 
lGvr [Ms)? S n0 + -jf-ipaW&H*)} (3-40) 
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are polynomials of degree 3 — e a in s. (Remember that <fi(s) is proportional to s — 4M 2 .) 
Consequently, 



W a (s) = w ^\ s ) K n(s) + polynomial + 0(p 

n=0 



8-2ea^ 



where K n (s) denote the dispersion integrals 



KM) 



dec k^iyOc^j 



7T X X — S 



(3.41) 



(3.42) 



The formula ( |3.41| ) contains and extends the one-loop result ( |3.13| ): By definition, 

K (s) = J(s) , (3.43) 

and the leading 0(p 4 ) part in w^(s) which contributes to Eq. ( |3.13j ) is manifest in the expres- 
sion (137401 ). 

As for the polynomial part of Eq. ( |3.41| ), its degree in s can be at most 3 — e a : The 
coefficient of any power of s higher than 3 — e a should not blow up in the limit M n — > 0, s fixed, 
and consequently, any such power can be absorbed into the neglected contribution 0(p 8_2£a ). 

The functions W a (s) then satisfy the bounds (|3.7D, since the functions K n {s) grow at 
most as (Ins) 2 . This means that Eq. ( |3.41| ) determines the amplitude A(s|t, u) up to a general 
crossing symmetric polynomial 5 A(s\t, u) of at most cubic order. The most general polynomial 
of this type contains six arbitrary parameters and it can be written as 



5A(s\t, u) 



+ 



+ 



(t - 2M 2 ) 2 + (u- 2M, 



2\2 



Ml 5(3 ( 4 „ 

7T 7T 

4(« " 2M 2 ) 3 + ^[{t- 2M 2 ) 3 + (u- IMlf 

"it "it 



(3.44) 



This formula represents the polynomial contribution to A(s|t, u) beyond the one-loop order 
( 3.17 ). In the chiral limit, the six constants involved in Eq. ( |3 . 44 ) dominantly behave as 



6a, 5(3 = 0(p 4 ) , SX li2 = 0(p 2 ) , A 3 , 4 = 0(1) 



(3.45) 



Furthermore 5a, 5(3, 5X± and 5X 2 can be absorbed into the parameters a, (3, Ai and A2, re- 
spectively, which characterise the one loop result fl3.17j ). This redefinition simply extends the 
expansion ( |3.19| ) in powers of quark masses and consequently the parameters a, (3 and Ai^ so 
redefined will be referred to by their original names. Notice that in the expressions ( |3.41|) for 
the functions W a (s), the redefinitions a — > a — 5a, (3 — > (3 — 5(3 and Ai 2 — > Ai j2 — 5\\^ induce 
terms of higher chiral order which can be reabsorbed into the neglected 0(p 8_2ea ) remainders. 
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This implies that the whole tttt scattering amplitude up to and including two loops can be 
expressed in terms of only six parameters, and the result reads: 



A(s\t, u) 



s - -Ml + a 



3 

2Ml? 



Ml 



3FI 



+ 



A2 
Ft 



(t - 2M^Y + (u- 2M, 



2\2 



+ ^{s-2Mlf + ^[{t-2Mlf + {u-2Mlf 
+ K(s\t, u) + 0(p 8 /A 8 H ) , 



(3.46) 



where 



+ 



+ 



n=0 



1 r 



K(s\t,u) = 32nJ2\^^\s)-w^(s) K n {s) 



.(")/ 



i[4 n) W + 3( S -«M n) (t)] K n (t) 



2 
1 r 



wf\u) +3(s-t)< j (u) K n (u) 



(3.47) 



The one-loop contribution Jr a ^\(s\t, u), cf. Eq. (|3.18|) , is hidden in the n = part of 
the sum (|3.47|) : It arises from the first terms, proportional to i^ 4 , in (|3.40|) . The terms 
proportional to F~ 6 in the polynomials w^\s) represent the 0(p 6 ) (and 0(p 7 )) contributions: 
The parts linear in Ai and A2 correspond to single loop graphs with one 0(p 4 ) vertex, while the 
remaining parts, which are cubic in a and (3, give rise to the genuine two-loop contributions. 

If one wishes to identify a given chiral order in the general formula ( |3.46| ), the parameters 
a, P, Aj, i = 1,2, 3, 4, have to be expanded in powers of quark masses. The only difference with 
the previous expansion ( |3.19|) , is that now more terms have to be included, until one reaches 
the first neglected order 0(p 8 ). Hence 



a 



n=0 



/?=£ P {n) , A 1>2 = £ A^ 



(n) 



2 ! 



A 



3.1 



n=0 



n=0 



2^ A 3,4 , 
n=0 



(3.48) 



where = 0(p n ) and likewise for the other constants. The particular case of the standard 
%PT is again characterised by the absence of odd chiral orders. 

Finally, we turn to the dispersion integrals K n (s) fl3.42| ) which occur in the two loop 
formula ( [3.46 ), ( [3 .47 ). It is remarkable that all these integrals are elementary and that they 
can be expressed in terms of powers of the standard one- loop function J(s) ( 13.141) : 



22 



K (s) = J(s) 

s - 4M 2 - 1 

KM = _^tf lW -_ (349) 

*'W - T^7^{7^[^>)-f^ 2 [^>)- 2 ]}-S 

ft W = T^7^{ 16 ^>'- 2 + 8 ^W + ^^ + T}-3i + T52- 

The proof of these formulae and a discussion of some properties of the functions K(s) can be 
found in Appendix C. 



4 Phenomenological Aspects 



The preceding discussion has led us to the expression ( [3.460 for the nn amplitude at two- 
loop order, which depends on the six parameters a, (3 and A«, i = 1,2,3,4. A priori, neither 
the one-loop formula ( 3.17|) nor the two- loop expression ( |3.46|) , by itself, contains information 
about these parameters. On the other hand, Eqs. (|3.23|) , ( |3.24|) or (|3.25|) relate them to the 
low-energy constants, such as mBo, m 2 Ao, m£, L\, L2,..., which appear in the chiral expansion 
of the effective action Z cK . To lowest chiral orders, these same low-energy constants are also 
present in the expansion of other observables, such as the masses of the pseudoscalars, F n /Fx, 
the scalar form-factor of the pion, the K14 form-factors, etc... This kind of information, which 
reflects the Ward identities obeyed by the QCD correlation functions, is rather crucial and 
has always played a central role in applications of chiral perturbation theory. In particular, it 
allows one to pin down approximate values of the parameters (3, Ai and A2 from independent 
experiments, and it provides a clear interpretation of the leading order value of a in terms of 
the mechanism of spontaneous breakdown of chiral symmetry in QCD, by relating this value to 
the quark mass ratio r = m s /m, cf. Eq. ( |3.22j ). It is also clear, however, that the higher orders 
in the expansions of a, (3 and the Aj's in powers of quark masses involve a rapidly increasing 
number of new low-energy constants from Z eS . The task of estimating them from independent 
measurements of different observables, and of re-evaluating the values of the previous low- 
energy constants taking into account the higher-order corrections, seems rather out of reach at 
the moment. 

Fortunately, it turns out to be possible to pin down four of the parameters (Ai, A 2 , A 3 
and A 4 ), which enter the two-loop expression (|3.46|) of the 7T7t scattering amplitude, using sum 
rules and available titt data at medium energies. The derivation of these sum rules and their 
evaluations are sketched in the following subsection (a more detailed account will be postponed 
to a forthcoming publication). Next, we consider a few applications. In particular, we discuss 
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how the low-energy data constrain the remaining two parameters a and (3. We remind the reader 
that a is the key parameter for deciding whether the chiral symmetry is broken according to 
the standard scenario, with a large (qq) condensate, or not. 



4.1 Sum Rule Evaluation of the A^'s 

Consider the s-channel isospin 1=0,1,2 amplitudes F I (s,t) defined as (u = AM 2 — t — s) 

F°(s,t) = — - {3A(s\t, u)+A(t\s, u)+A(u\s,t)} 

F\s,t) = -i- {A(t\s, u) - A(u\s, t)} (4.1) 

62,11 

F 2 (s,t) = -L{A(t\s, u)+A(u\s, t)} . 

Because of the Froissart bound, they satisfy twice-subtracted, fixed-t, dispersion relations. 
These form the starting point for the derivation of the sum rules for the constants \. The 
next, and somewhat tedious, step consists in transforming this set of fixed-t dispersion rela- 
tions into a form which allows direct comparison with the two-loop chiral expression. This is 
performed by first identifying and rejecting the contributions which are of chiral order 0(p 8 ) 
or more. In order to do so, one has to split the integration ranges of the dispersive integrals 
into two regions, a low-energy region [4M%, E 2 ] and a high energy region [E 2 , oo]. In the first 
region, one may drop the contributions from the partial waves with I > 2, while in the latter 
region one may expand in powers of s, t, u divided by E 2 and stop at the third power. The 
borderline energy E corresponds, roughly speaking, to the point where the chiral expansion 
starts to break down, i.e. E is in the range 500-600 MeV. Secondly, one has to impose crossing 
symmetry. In practice, one proceeds in close analogy with subsection III-C of Ref. ||. There, 
starting from thrice-subtracted dispersion relations, it was shown how, up to corrections of 
chiral order 0(p 8 ), crossing symmetry entirely determines the form of the amplitude except 
for six arbitrary parameters. If, instead, one uses twice-subtracted dispersion relations one 
finds that two parameters remain undetermined (a and j3), while the four others (Ai,...A4) obey 
well- convergent sum rules 4 . Details will be provided elsewhere. Here, we shall only give the 
final result: 

^ = ~ {4 2) - 4 2) ) + 2M 2 (Zf - if) - 3X« - mil? + ^(h - 4h 2 ) 



^ = I(zf_4*>) + zf> + ^(-9/^-/4) 



4 As a matter of fact, if one assumes conventional Regge asymptotics for the t-channel 1=1 amplitude, one 
can derive a slowly convergent sum rule for (3 as well, which is equivalent to the one derived long ago by Olsson 
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4 = V - V + — (K - K) 

F 6 2 2 9 2 



(4.2) 



In these formulae, the quantities are given in terms of integrals over the low-energy range 
[AM 2 , E 2 ], while h n and h' n involve integrals from E 2 to infinity (detailed formulae will be given 
below). Since experimental data on tttt phase shifts exist down to roughly 500 MeV, we can 
calculate the latter integrals essentially using experimental input. In the low-energy integrals 
we shall use the two-loop chiral expression for the phase shifts. Stability of the result with 
respect to variations of the energy E within the range indicated above measures the extent to 
which the phases do match to the experimental ones at this energy. Note that since the 
%PT phases themselves depend on the parameters Aj, the set of relations Q4.2|) must actually 
be solved in a self-consistent way. 

Let us now give the explicit expressions for the various entries in Eqs. ( |4.2|) . First, we 
define the Zl fc )'s as 



32vr d 
k\ ds k 



71 



dx ImWn 



AM 2 X 



2-e a 



- E^ n) (^n(* 



X 



(4.3) 



n=0 



where the polynomials w a n \s) and the functions K n {s) are those which occur in the expression 
of the two-loop amplitude of the preceding section. The quantities h a and h' a are defined in 
terms of the following high-energy integrals 



, , , -1 f» , 2x + t - AM 2 
h>o(t) = — / dx— — 7— r-rrlm 

v 1 tt Je* x 2 (x + t-AM 2 ) 2 



±F°(x,t)+F\x,t) + ^F 2 (x,t) 



(4.4) 



and 



-1 f°° , 2x + t-4M 2 



— 1 r 00 
h 2 (t) = — / dx 

71 JE2 



x 2 (x + t-AM 2 ) 2 



Im 



l -F\x,t)- l -F\x,t) + l -F 2 {x,t) 



(4.5) 



Then, h a and h' a are given simply by 

h a = h a (Q) , 



0,2 . 



(4.6) 



Finally, hi is given by 



1 poo 

hi = — 



dx 



tt Je 2 x(x - AM 2 ) 



-Im 



V(*,0) + V(x,0)- V(x,0) 
6 2b 



(4.7) 



Notice that the three integrals (|4.4j) , (|4.5| ) and (|4.7|) involve the t-c/ianne/isospin b=0, 2 and 1 
combinations, respectively. In addition to the constraints it imposes on the low-energy parts, 
crossing symmetry also implies the following relation between these high-energy integrals: 



9/ii = 2h'r. — 5/io • 



(4.8) 
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In practice, the most reliable source of information on the tttt phase shifts remains the old 
high-statistics CERN- Munich experiment |34| (see e.g. the reviews by Ochs [|35| and by Morgan 
and Pennington [3E] for critical discussions of more recent experiments). In order to estimate 
the errors in the evaluation of the high-energy integrals, we have varied the experimental data 
within their error bars and we have also compared the results obtained using different analyses 
of the production data of Ref. [33J]. The errors in the low-energy integrals, which reflect the 
presence of 0(p 8 ) contributions, were estimated from the sensitivity to variations in the value 
of the matching energy E and by trying different ansatze for Im/ a (which all differ only by 
0(p 8 ) contributions) in equations (fO|). Furthermore, a small dependence of the result upon 
the remaining two parameters a and j3 was absorbed in the error. More details of this analysis 
will be given elsewhere. We end up with the following estimates for the values of Ai and A2: 



Ax 



-5.3 ±2.5) x 10" 



A 2 = (9.7 ± 1.0) x 10" 



(4.9) 



These numbers are in agreement with the result of the analysis of Bijnens et al. [3~I|, which 
incorporates constraints from sum rules for D-wave threshold parameters as well as from K14 



decays. An independent estimate was recently made by Pennington and Portoles ||37|| , who 
derived simple, but approximate, formulae by setting the isospin 1=2 contributions in both 
direct and crossed channels exactly equal to zero. From their analysis it follows that 
Ai = (-4.3 ± 2.3) x 10" 3 and A 2 = (8.8 ± 1.1) x lO" 3 for i^=92.4 MeV. Our result for A 2 lies 
between those of these two references. Finally, for the remaining two parameters, we obtain 



A 3 = (2.9 ±0.9) x 10" 



A, 



-1.4 ±0.2) x 10" 4 . 



(4.10) 



4.2 Low-energy Phase Shifts 

It is well known that the chiral expansions of the nir partial- wave amplitudes do not satisfy 
unitarity exactly, but rather in a perturbative sense. The deviations from unitarity provide an 
estimate for the energy range where the chiral expansion, at a given order, can be trusted. We 
show in Fig. 1 the Argand diagram for the amplitude /q computed at one-loop and at two-loop 
order. The figure shows that there is a considerable improvement with respect to unitarity at 
low energies in going from the one- loop result to the two- loop result. On the other hand, it 
also suggests that one should expect sizeable 0(p 8 ) corrections starting at 450 MeV. 

In spite of the fact that the unitarity relation is not exactly satisfied, it is possible to 
define a chiral expansion for the phase shifts (see, for instance, the discussion in [0). At 
lowest order, and corresponding to the amplitude ( |3.9|) , the partial waves for / < 2 are defined 
as (a = 0,1,2) 

/ o _ 4 M 2 

Sa = \ + 0( P 4 ) . (4.11) 
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The one-loop form ( |3.17| ) of the amplitude extends the above definition of the phases to 

5a - 



! s -AM 2 



[<Pa + 1>a] + 0(p b ) 



(4.12) 



whereas the two-loop expression (|3.46| ) leads to 

8a = \ 



4M 2 



Re f a + -((f a 



0(p 8 



(4.13) 



where in this last formula Re/ a must be evaluated from the full two- loop amplitude ( |3.46| ). 
The higher partial waves receive contributions starting only at the one-loop approximation, 
otherwise the formulae are exactly the same. The result for the phases as computed from the 
successive approximations ( 4.11 ), ( 4.12 ) and ( [4.13 ) are shown in Figs. 2a and 2b. For simplicity, 
we have used the same values for a and (3 (a = 2, (3 = 1.08 ) in all three cases (strictly speaking, 
of course, one should use the expansions ( ET48D of a and (3 at the given order). Surprisingly, 
for 1 = 1 = (Fig. 2a) the one-loop and two-loop phases remain fairly close at energies where 
one would expect (say from the Argand diagram) that 0(p 8 ) terms should become important 
(approximately above 500 MeV). This is somewhat fortuitous, and does not happen for the 
isospin 1=2 wave or the P-wave. In the latter case, it is interesting to observe how the p 
resonance is gradually built up as the order of the perturbation increases. Finally, Fig. 4a 
shows the behaviour, at increasing orders of the chiral expansion, of the difference 8$ — 5\. 

We investigate next the sensitivity of the phase shifts on the values of the parameters 
a and (3. Results are shown in Figure 3, representing 5q, 5q and 6{ respectively, and Fig. 4b 
which shows the difference 5® — 5{ in the energy range accessible in decays. We have taken 
values of a ranging from 1 to 3. Figure 4b suggests that the data of Rosselet et al. 
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arc 



approximately compatible with such a range. The figures show that the sensitivity upon a is 
most significant at very low energies. An important issue for us is to assess the extent to which 
the availability of a new generation of data on Kja decays from e.g. the DA$NE facility will 
allow a really precise determination of these parameters a and (3. For this purpose we show 
in Fig. 5 how the errors on Ai,...,A4 affect the prediction for 5$ — 8\. The figure shows that 
the situation where a = 1 and the one where a = 2 can be distinguished, provided sufficiently 
precise data are available in the energy range below E = 340 MeV. Let us now examine in 
a more quantitative way how the presently available data constrain the values of a and f3, 
using the two- loop form of the amplitude and our estimates of the values of the parameters Aj. 
Since there are only five data points provided in Ref. 



39 , it is useful to look for an additional 



constraint between a and (3. Such a relation can indeed be obtained from the "Morgan-Shaw 
Universal Curve" [ |D| . The latter represents a correlation between the 1 = and 1 = 2 S-wave 
scattering lengths, obtained by fitting the data in the region between 600 MeV and 1 GeV (i.e. 
outside of the range of applicability of the two-loop chiral expansion) using the Roy equations 
[3] ). We reproduce here the form quoted by Petersen in Ref. |24|| : 



(see fll 



2a° - 5a 2 = 0.692 ± 0.027 + 0.9 (ajj - 0.3) + 1.2 (a° - 0.30) 



(4.14) 
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Upon using the explicit expressions for the scattering lengths (see Appendix D), this equation 
transcribes into the desired relation between a and ft, which is plotted in Fig. 6. We have 
performed a chi-square minimisation using this relation in addition to the data of Rosselet et 
al. [0 and to the results ( PD , flOPP . We find 



a = 2.16 ±0.86, ft = 1.074 ± 0.053 , (4.15) 

with a x 2 equal to 2.01 for four degrees of freedom. The threshold parameters will be discussed 
in detail below. Let us give here the corresponding value for a®, 

a° Q = 0.263 ±0.052 . (4.16) 



Interestingly enough, this result is in agreement with the one quoted in Ref. |p4}| , which is 
based on a fit using the Roy equations. 

A physically interesting quantity which is related to the tctt phase shifts is the phase of 
the CP violation parameter e' 

0( e ') = 90 + tfjj - 5° (degrees) , (4.17) 

where the phases are to be evaluated at the energy E = M K o = 497.7 MeV. The preceding 
discussion shows that it makes sense to use the two-loop approximation at this energy and, 
indeed, Gasser and Meifiner have discussed 0(e') at 0(p 4 ) in the standard chiral expansion and 
find = (45 ± 6)°. Using the values of a and ft (|4.15|) derived from the fit discussed above, we 
obtain 

<p= (43.5 ±2 ±6)° , (4.18) 



where we have shown separately the error coming from the uncertainties in the fit values ( 4.15|) 



of a and ft (±2°), and the error due to the uncertainties in our determinations ( |4.9| ) and (|4.10|) 
of the Aj's (±6°). At this energy, a significant gain in precision is achieved by going through a 
numerical solution of the Roy equations, rather than by using the chiral expansion. This has 
been done by Ochs , who found a result very similar to ours but with a much smaller error: 
(j) = (46 ± 3)° . 

4.3 Threshold Parameters 



From the two-loop amplitude ( 3.46|) , it is straightforward to obtain explicit expressions 



for the scattering lengths a\ and the slope parameters bj. For the reader's convenience, we have 
gathered, in Appendix D, a selection of formulae relevant for the discussion that follows. 

Table 1 shows the numerical values of the threshold parameters at two-loop precision for 
different values of a and ft, with the Aj's fixed at their central values as given by Eqs. (|4.9|) 
and (|4.10|) above. The second line of the table gives the two-loop threshold parameters for 
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a 




0° 

u 




-10a 2 , 


-1068 


1.04 


1.08 


0.210 


0.26 


0.43 


0.78 


2 


1.127 


0.264 


0.27 


0.32 


0.83 


1.071 


0.255 


0.25 


0.29 


0.78 


1.016 


0.246 


0.24 


0.26 


0.74 


2.16 


1.074 


0.263 


0.25 


0.27 


0.79 


2.5 


1.085 


0.283 


0.25 


0.22 


0.80 


3.0 


1.102 


0.311 


0.25 


0.16 


0.82 


3.5 


1.122 


0.342 


0.25 


0.10 


0.85 


Error bars 


±0.006 


±0.02 


±0.02 


±0.06 


Experiment |24 


0.26±0.05 


0.25±0.03 


0.28±0.12 


0.82±0.08 



a 




10a} 


10 2 b\ 


10 2 a° 


10 3 a 2 


io 4 4 


1.04 


1.08 


0.37 


0.60 


0.17 


0.08 


0.56 


2 


1.127 


0.39 


0.55 


0.17 


0.13 


0.65 


1.071 


0.37 


0.55 


0.17 


0.14 


0.63 


1.016 


0.35 


0.55 


0.18 


0.14 


0.62 


2.16 


1.074 


0.37 


0.54 


0.17 


0.15 


0.65 


2.5 


1.085 


0.37 


0.53 


0.17 


0.17 


0.68 


3 


1.102 


0.38 


0.50 


0.17 


0.20 


0.74 


3.5 


1.122 


0.38 


0.48 


0.17 


0.23 


0.80 


Error bars 


±0.01 


±0.15 


±0.05 


±0.28 


±0.18 


Experiment |24| 


0.38±0.02 




0.17±0.03 


0.13 ±0.30 


0.6 ±0.2 



Table 1: Threshold parameters for F„ = 92.4 MeV, in units of M n + = 139.57 MeV for 

different values of a and of (3. 
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the values of a and (3 corresponding to the 0(p 4 ) predictions ( |3.28| ) of the standard case f\ 
within the SU(3)l x S77(3)r analysis, while the fourth line corresponds to the central values 
( J4.15|) , obtained from the fit of the two-loop phases to the Rosselet data cum Universal Curve 
constraint. For the remaining entries, we have, for each value of a, chosen for (3 the central 
value given by the Universal Curve (|4.14j) , except for a = 2, where we have also displayed 



the results corresponding to the maximal and minimal values of (3 as allowed by ( f4.14| ). This 
variation of (3 at fixed a only affects the values of the S- and P-wave scattering lengths and 
the S-wave slopes, which are the only threshold parameters to depend on (3 already at leading 
order. The effect is strongest in Oq and is of the order of about ±10% in average. The S-wave 
scattering lengths a° and are clearly the most sensitive to variations in the value of a. The 
penultimate line of Table 1 gives the error bars (which are independent on both a and /3 at the 
level of numerical accuracy we consider here) on the threshold parameters. They only include 



the effects of varying the values of the Aj's within their ranges as given by Eqs. (4.9) and ( 4.10 



The last line gives the values, as quoted in |2J], obtained from Roy-equation analyses of the 
and production data. 

Figure 7 shows the S-wave scattering lengths as functions of a. The shaded area gives 
the uncertainties coming from the errors on the Aj's, as given by Eqs. (|4.9| ) and (|4.10Q , whereas 
the band delimited by the solid lines includes the error on the a-dependence of (3 inferred from 
the Universal Curve Q4.14 ). The present experimental value a° = 0.26 ± 0.05 allows a to lie 



between 1 and 3.4. In this whole range of a, corrections to a[j represent a 25% increase with 
respect to tree level, and the two- loop corrections are about 5% of the one- loop result, (see also 
the discussion at the end of Appendix D). The observed convergence rate leads one to expect 
that higher order effects should be within the error bars quoted in Table 1. In addition, we 
show (Fig. 7) the combination \a^ — a^\ which can be extracted from the measurement of the 



7r + 7r -atom lifetime planned at CERN flQfl . 



5 Testing the Strength of Quark Condensation 



In Section 4.2 it has been shown how the two-loop expression for the scattering amplitude 
can be exploited to extract values of the parameters a and j3 from available low-energy data. 
The lack of precision in the result, viz. a = 2.16 ± 0.86, merely reflects the large error bars 
in currently existing data, rather than the uncertainty attached to the parameters Ai, . . . , A4. 
More accurate data are awaited, and with them a rather precise measurement of a and (3 should 
become possible. We have already stressed the key role played by the parameter a in the 
measurement of the amount of quark condensation in the QCD vacuum: In the standard, large 



5 Of course, the two-loop values of a s t and (3 s t need not coincide with their one-loop evaluations (3.28). 
However, the changes are expected to be small, and we feel that the second line of Table 1 gives a realistic 
description of the standard two-loop predictions within the error bars quoted in the penultimate line. 
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condensate alternative of symmetry breaking, a should remain close to 1 (cf. Eqs. ( |3.28| ) and 
(|3.33|) ), whereas the case of a marginal quark condensation is characterised by a considerably 
higher value, a ~ 3 — 4. In both cases, ft should stay close to 1. 

In this section, we attempt a more quantitative interpretation of the low-energy constants 
a and ft, relating them to fundamental QCD parameters such as the quark mass ratio r = m s /m 
and the condensate parameter 2m_B - To leading order, this relation has been given before 
|| p. The complexity of this relationship rapidly increases with increasing chiral order to 
which a and ft are expanded. Here, we restrict our discussion to one- loop order 0(p 4 ), using 
the corresponding perturbative expressions Q2.171 ) and ( |2.18| ) for ft and a. Let us stress that this 



is a perfectly consistent procedure: It is conceivable that two-loop accuracy in the expansion 
of the amplitude is important for a precise measurement of a and ft using, in addition to other 
input, energy- dependent phase shift data above threshold. On the other hand, to make use of 
the results of such a measurement, one- loop accuracy in the expansion of a and ft in powers of 
quark masses could be sufficient to distinguish a value of r ~ 25 from, say, r ~ 10, or the case 
2mB /M% ~ 1 from 2mB /M% ~ 0. 

5.1 Dimensional Analysis and Order of Magnitude Estimates 

Before proceeding to a more refined analysis, it might be useful to recall some crude and 
simple order of magnitude estimates of various contributions to a and ft that are exhibited 
in Eqs. (2.17) and (2.18). A constant of C eS multiplying an invariant that is made up from 



n covariant derivatives and m powers of the quark mass is, for n + m > 2, expected to be 
of order F$ A^~ n_m , unless it is suppressed by the Zweig rule. Indeed, this is the expected 
order of magnitude of a (massless) QCD correlation function that consists of n vector and/or 
axial- vector and of m scalar and/or pseudoscalar current densities, at low external momenta, 
with the singularities arising from exchanges of Goldstone bosons subtracted. This estimate 
is obtained by saturating correlation functions that are smooth enough at short distances by 
the lowest massive bound states (p, ax, fo, • • •), which have masses of the order of Ah ~ 1 GeV. 
(F is a scale characterizing the coupling of mesons M to currents J: (0\J\M) ~ F A H .) This 



argument is a special case of the "naive dimensional analysis" of Ref. fiB|] , except that we leave 
open the question as to whether it may be extended to the "one-point function" (n — 0, m = 1) 
representing the quark condensate. The reason for this restriction is quark confinement, which 
prevents us from analysing the quark condensate by the same dispersion techniques as in the 
case of (n + m)-point functions for n + m > 2. 

Factoring out F 2 from all terms of C eS , cf. Eqs. (PD , (gig ), ( gig ) and (gig) , the £ (2) 
constants A and Zq are two-point functions divided by Fq, and consequently they are expected 
to be of the order of 1. (A more detailed estimate, using QCD sum rules, gives A ~ 3 ± 2.) 
Similarly, the constants £, pi, p 2 are expected to be of the order of l/A H , whereas the 
constants Ai, A 2 , A 3 , Bi, B 2 , Ex, E 2 , E 3 should be of the order 1/A 2 ^. The constants Zq, £, p±, 
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p 5 and pt violate the Zweig rule, and they are therefore expected to be suppressed compared 
to the corresponding dimensional estimates. 

It should be stressed that for n + m > 2 these estimates are independent of any chiral 
counting or of the importance of the quark condensate. On the other hand, attempts to 
estimate the value of the u and d running quark masses f], especially rh = {m u + m,d)/2, depend 
in one way or another on the presumed mechanisms of chiral symmetry breakdown in QCD. 
A model-independent evaluation of corresponding QCD sum rules is at present problematic, 
due to the complete absence of experimental information on the size of the spectral function of 
the divergence of the axial current beyond the one-pion contribution. The existing evaluations 
of m usually normalise [15] the unknown spectral function by its low-energy behaviour as 
predicted by the standard xPT. However, in generalised xPT the latter prediction is modified 
by an enhancement factor which can be as large as 13.5, depending on the actual importance of 



the quark condensate [46], |17| . The issue can and should be settled experimentally, by measuring 
the magnitude of the divergence of the axial current in hadronic r decays [f5| . Fortunately, the 
existing determinations of the quark mass difference m s — m u (see e.g. the recent analyses in 
|}47| , f48| , |49| , from which earlier references can be traced back) are on less speculative grounds. 
They are based on sum rules which involve the two-point function of the divergence of the 
vector current sj^u; the corresponding spectral function can be normalised using experimental 
data in a rather model-independent way. For definiteness, we shall adopt the corresponding 



result of Ref. ||47 |, viz. m s — m u = (184±32) MeV, which is compatible, within the error bars, 



with the values obtained by the authors of Refs. IE8J and H3|. For the sake of our order of 



magnitude estimates we shall thus use for m a value which depends on the ratio r = m s /m in 
the following way 



m = - ^MeV, 5.1 

r — 1 

keeping in mind that in generalised ;\RT, r is not fixed from the outset. Even for r as small as 
r ~ 8, the corresponding value of m would be 

m ~ (26 ± 4.6) MeV , (5.2) 



which would keep the expansion parameter fn/An reasonably small, although somewhat larger 
than in the standard %PT. Given the above value for m s — m^, Eq. ( |5.2| ) represents the upper 
bound for m in G^PT. 

We are now in a position to estimate orders of magnitude of various terms contributing to 
a, dividing Eq. ( |2.18j ) by M%. For definiteness, we set A H ~ M p and we first consider the low B 
alternative, taking, e.g., r ~ 10. In this case, both terms TmB^jM'l and 16m 2 Ao/M% should 

6 All QCD running quantities will be understood to be normalised at the scale 1 Gev ~ Ah- The combinations 
of constants times powers of quark masses appearing in the expansions of a and (3 are QCD renormalisation 
group invariant. 
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be of order unity. The £® contributions 4m£ and 81m 3 p 1 /M% should both reach the ten- 
percent level, whereas the £(2,2) terms m 2 Bi, m 2 A 3 and the £(o,4) terms, including their large 
coefficients - viz., 256m 4 Ei/M 2 - should hardly reach the one-percent level. Chiral logarithms 
and the Zweig rule violating, (explicitly) r-dependent, terms will be discussed shortly. 

In the standard large Bq alternative (r ~ 26), the whole expression ( |2.18| ) is dominated by 
the Gell-Mann-Oakes-Renner (GOR) ratio 2mB /M 2 , which remains close to 1. All remaining 
terms should contribute at most at the one-percent level, reflecting the smaller value of m ~ 6 
MeV. 



5.2 Constraints from Pseudoscalar Meson Masses 
and Decay Constants 

More information on the low-energy constants contained in and can be obtained 
from the expansion of the Goldstone boson masses and decay constants F„ and Fk- These 
expansions, up to and including order 0(p 4 ) contributions in generalised x?T, are collected 
in Appendix A. To leading, 0(p 2 ) order, m 2 Ao and rr^Z^ can be expressed as functions of 
r = m s /m and of pseudoscalar meson masses: 



4m 2 A r 2 -r _ 

M 2 ~ r 2 - 1 ~ e[r) 

(5.3) 

8m% p Agmq 



M 2 (r ~ If 

where (using M v = 135 MeV, M K = 495.7 MeV, M v = 547.5 MeV 



r 2 = 24-1 = 26 
Ml 



Agmo = (3M 2 + M 2 -4M^)/M 2 = -3.6 



The ratio 2mB /M 2 depends, in addition, on the Zweig rule violating constant Zq\ 

2mB 



Ml 

where 



(5.4) 



e(r)[l + (r + 2)C] , (5.5) 



7 s T 

^0 -^8 



is expected to be small. For r = r 2 + 0(m quark /A H ), the ratio ( |5.5|) is close to 1, whereas the 
constants (|5.3|) are relegated to higher orders, provided, of course, that Agmo = 0(m quar k/Ajy). 
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This is the standard scenario, in which the quark condensate is large enough to dominate the 
expansion of Goldstone boson masses as well as the expansion ( |2.18|) of the parameter a. As r 
decreases, the ratio ( |5.5| ) becomes smaller than 1 and (neglecting £) it vanishes for 

r = n = 2—^ - 1 = 6.3 . 5.7 
Notice that r cannot be smaller than this critical value, since vacuum stability does not allow 



the ratio ( |5.5| ) to be negative. For r < r 2 , the ratio e(r) ( |5.3|) gradually rises towards 1, in 
agreement with the dimensional analysis mentioned above: Assuming the quark mass m to be 
given by Eq. ( |5.1| ), and taking, for instance, r = 10, Eq. ( |5.3| ) implies e ~ 0.3 and Aq ~ 3.4. 
In generalised xPT, the leading order contribution to a is obtained by summing up the first 
three terms in the expression ( [2.18|) : 



a<« = l + 6^-4(l + 2C) , (5.8) 
r z — 1 

whereas 

/3 {0) = 1 . (5.9) 

Hence, the quark mass ratio r, instead of being determined by pseudoscalar meson masses, may 
be used as a measure of the amount of quark condensation. It determines the leading order 
value of the parameter a, which lies in the range 1 < < 4. 

Let us now consider the constants of characterizing the next to the leading order 
0(p 3 ). To that order, the constant m£ can be inferred from the splitting of the decay constants 
Fk and F n , cf. Appendix A (we take the determination Fx/ 'F n = 1.22 ± 0.01 from Ref. pUJ), 



1 (Fl \ 0.488 + 0.024 , r 1n . 

m « = —i (pf - V = —i ■ (5 - 10) 

This value agrees with the dimensional analysis estimate within a factor of 2. The next to the 
leading order contribution to (3 is then 

(5.11) 




The 0(p 3 ) corrections affect the leading order formulae ( |5.3| ), ( |5.5| ) and ( |5.8| ). The most impor- 
tant effect arises precisely from the splitting F K ^ F n . It amounts to replacing e(r) everywhere 
by e*(r), where 



r* 2 -r » n (F K M K 



r 2 - 1 ' z V F w M n 



6*(r) = 2^— , r 2 * = 2 -^M -1. (5.12) 



2 



The relative importance of this correction is greatest for r ~ r 2 (r 2 = 26, whereas r 2 = 39). In 
the low B alternative, this £(2,1) correction shows up mainly as a slight increase in the critical 
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value of r for which the quark condensate vanishes: This critical values moves from r = t\ ~ 6.3 
towards r = r*, with 

Finally, we discuss the £(0,3) parameters pi, especially p\ and p 2 which are not suppressed 
by the Zweig rule. The magnitude of one of them - the parameter p 2 - can be inferred from 
isospin breaking effects induced by rrid 7^ m u . It has been pointed out in || that (in the absence 
of electromagnetism) there exists a particular linear combination of the K + , K° and ir + masses, 
namely 

A(r, R) = M 2 K - Ml - RAMI + « - "A^K) ^ , (5-14) 

2 r + 1 



where 



nrn — nrn 

R=— , AM 2 K = (M 2 R0 - M 2 K+ ) QCD , (5.15) 

m d -m u 



which has remarkable properties: Neglecting terms quadratic in — m u , all contributions to 
A(r, R) of the type m^^Bo, Wq Uark and even In m quar k cancel out. Consequently, in 

the standard ;\TT, A(r, R)\ md=mu = O (m^ uark ) . In generalised xPT, A(r, R) also vanishes at 
leading order; however, it receives 0(p 3 ) contributions from C^ 3 \ In both cases, 



A(r, i?)| Wd= ^ = I - 1 I (M£ - ^-M^j + (r - l)\r + l)m> 2 + O(m^ uarkJ B , m* uark ) . 

(5.16) 

Let us write, as usual, 

AMI = (M 2 K o - M| + )ex P + l(Ml + - M 2 ), (5.17) 

where the coefficient 7 — 1 measures the violation of Dashen's theorem. We shall take 

7 = 1.8 ±0.4, (5.18) 



based on the recent estimates [51]. Taking for the ratio R the value 52 



R = 43.7 ±2.7 (5.19) 
inferred from baryon mass splittings and from u — p interference, one obtains 

14.4 ±2.1. (5.20) 



RAM 2 K 



Ml 



Upon inserting this value into Eq. ( |5.16|) , one finds that the central value of p 2 is largely 



dominated by the uncertainty: For r < 15, p 2 is compatible with zero; more precisely, 

-0.6 ±2.1 

Ml (r- l) 2 (r + 1) ' K ' } 
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In comparison, for r = 25, the above numerator would be —2.9 ±2.1. It is interesting to notice 
that the dimensional analysis of Section 5.1 gives (the error bars come from the uncertainty on 

0.4 ±0.2 



m 3 pi 



Ml 



DA 



D 3 



1,2 



(5.22) 



The constants p\ and pi also describe the leading SU (3)v-breaking effects in two-point functions 
of the scalar and pseudoscalar densities, which might lead to additional information. We have 
however not performed a complete analysis along these directions so far. In the sequel, we shall 
use the estimate ( |5.22j ), since it is compatible with the more direct determination from A(r, R) 
in the case of p 2 . 



5.3 a and (3 as Functions of the Quark Mass Ratio r = m s /m 

Up to now, nothing has been said about chiral logarithms. In generalised xPT, they start 
to contribute at the order 0(p 4 ), i.e. in the expansions of a and of /3, they first appear in the 
next-to-next-to leading order corrections and (3^ 2 \ As already pointed out, one particular 
feature of generalised x?T is that 0(p 4 ) loops renormalise the 0(p 2 ) constants Aq and Zq 
by higher order terms proportional to Bq (remember that Bq counts as O(p)). Similarly, the 
constants £, £ and of are also subject to 0(p 4 ) renormalisations, proportional to Bq. The 
constants of O 2 ^ and of £y\ which appear in the expressions Q2.18|) and ( |2.17|) for a r and for (3 r 



(or in the formulae of Appendix A), are the renormalised constants at scale p, with the 0(p 4 ) 
So-dependent counterterms exhibited in Eq. (|2.11|) included. For a quantitative discussion 



of 0(p 4 ) effects, it is useful to explicitly split off these //-dependent parts from the constants 
that are dominantly of order 0(p 2 ) or 0(p 3 ). This separation is needed in particular when 
dealing with the Zweig rule violating constants Zq, £, p 4 , p 5 and pj\ While it makes perfect 
sense to declare these constants to be small at their respective 0(p 2 ) and 0(p 3 ) tree levels, 
they will, however, be generated at the 0(p 4 ) accuracy by Zweig rule violating chiral loops. In 
what follows, we shall assume that there exists a scale A (not related to the renormalisation 
scale p) at which these constants vanish. Upon using the coefficients Tx, Eq. ( |2.15|) , of the 
corresponding beta functions 7 , this assumption amounts to writing 

Am 2 Z3 11M 2 - 2 



o x2ln ^ 



Ml 288tt 2 F2 A 2 

~~ K ' p 2 
- 3 p 4 mi x( ny_ 2z]l X 



Ml 192tt 2 F 2 V 3 / A 2 



7 For Zq and for £, they can be inferred from the renormalisations of the constants Lq and L4, respectively, 
as given in Ref. ||, whereas for the pt's, see Ref. 28 . 
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Ml 647r 2 F% V 3 J A 2 

m 3 p 7 Ml , /i 2 
= tz n — 

Ml 576vr 2 F 2 A 2 ' 

where we have introduced the notation 

_ 2mF 2 B 4m^VA) _ 4m 8 JgZf 

F 2 M 2 > ^ F 2 M 2 > F 2 M 2 ■ 

The vanishing at a common scale A of all the and Zweig-rule violating constants is 
by no means obvious. Here it is assumed mainly for reasons of simplicity. An experimental 
determination of the sizes of Zq(A) and of £(A) is in principle possible from a simultaneous 
analysis of K14 form factors and of the low-energy 77 — > 7r°7r° cross-section, when more precise 
data on these quantities become available 0. 

It is now rather natural to use the scale A in order to separate the induced 0(p 4 ) compo- 
nents of the constants Aq(ji) and £(p). In terms of the notation introduced in Eq. Q5.24Q , one 
has 

4m 2 F 2 A (^ 5M 2 2 p 2 

FlMl V 96ir62F? X n A 2 

7T 7T 7T 

(5.25) 

3M 2 u 2 
2m£ = ^(AJ-g^sln^. 

Similarly the constants pi(p) and P2(p) of £^ 3 ^ which are not suppressed by the Zweig rule 
may be decomposed as 



F*M* r 64vr 2 F 2 V 6 J A 2 

(5.26) 

m 3 F 2 p 2 ^ Ml fxy \ p 2 



F 2 M 2 rz 64vr 2 F 2 V 6 /A 

where 



m 3 F 2 p,(A) 



F 2 M 2 



1, 2 . (5.27) 



We are now in a position to reexpress various contributions to a and (3 in terms of 
observables, including 0(p 4 ) accuracy. From the expressions of F 2 /Fq and F^/Fq given in 
Appendix A, one obtains 



r(2,2) r(2,2) 
°F,K ~ °F,ir 



3M 2 1 



+ x ln (^ 2 M 2 ) + - (-5p. + 2p^ + 3 N ) + 0(p 3 ), (5.28) 

6ln A b t r — 1 



8 An attempt to estimate the possible effects of non-vanishing Zq and £ has been made in a preliminary 
version of this work |26|. See in particular Fig. 1 in this reference. 
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where 8^'% — 8p^ collects all 0(p 4 ) tree contributions arising from £(2,2)- Of course, the entire 
expression ( |5.28| ) is //-independent. All results collected so far in this Section can now be used 
to reexpress the full parameter a, i.e. a r (/t) as given by Eq. ( |2.18| ), plus the logarithmic 
unitarity corrections displayed in Eq. ( |3.23|) , as follows: 



4 ( F 2 \ 

a = x + Ay - — - -1-1+81^ + ^ + A a (/i) + u a (x, y, r) + 0(p 3 ). (5.29) 



Here, A a (/i), which is of order 0(p 2 ), collects all the tree contributions arising from £(22) and 



£( 0j4 ). It can be read off from Eqs. (|2.18j) , (|5.28| ) and the formulae collected in Appendix A; 
however, its explicit expression would be of little use here. u) a (x, y, r) involves all the logarithms 
arising both from tadpoles and from the unitarity corrections: 



u a (x,y,r) 



M 2 rl 

-(1945 + 352r + 31r 2 )xy 
.3 



192tt 2 F 2 



+ (10 - Ar 2 )xz + -(148 + llr)x 2 - 24x 
o 



ln(/i 2 /A 2 ^ 



+ 

1 



10 

r — 1 
18 



- 3-292/ 



2jJL K 



r — 1 



+ 1 + (r + 7)y 



IN 

Ml 



r — 1 



+ 1 + I5y - 4(r - A)z 



32n 2 F? 



;i + 22j / + 33j/ 2 )[ln(M 2 //i 2 ) + l] 



Ml 



8 w^ 1 + ( r + Vvft 1 + 3 ( r + l )v\[HM z K /n + 1] 



Ml 



32n 2 F 2 



1 2 



1 + (2r + l)y 



r — 1 



-A 



GMO 



[ln(M 2 //i 2 ) + 1] 



(5.30) 



Similarly, the full parameter /3, i.e. (3 r of Eq. ( |2.17| ), plus unitarity corrections given in Eq. 
( gg| ), reads 



r - 1 I Fl 



(5.31) 



Here again, Aree(/-0; which is of order 0(p 2 ), collects all the tree contributions from £(2,2), 
whereas /3 loc , p (/x) contains the logarithms: 



Ml 



x ln(/x 2 /A 2 ^ 



+ 



16tt 2 F 2 
1 



r — 1 

M 2 



-5/v + 2/i X + 3/^,) 



16tt 2 F 2 



=-(2 + 5 2 /)[ln(M 2 //i 2 ) + l 



M. 



32n 2 F 2 



^[l + (l+r)y][HM 2 K /^) + l]. 



(5.32) 
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The parameters x and y ( |5.24| ) satisfy the following two equations, which may be obtained by 
multiplying the expansion formulae for FlM 2 and F^-M^ (see App. A) by the factor Fq/F 2 M 2 : 



x + y = 1 - 9pi - p 2 + A 7r (/i) +u w (x,y,mr) 



(5.33) 



r + 1 

£ H — - — y 



i(tsf)'- 3(r2+r+1) ' 5l - (r2 - r+1)ft 



2 r 

+ A^(/i) + u K (x,y,r) . 

As before, the /^-dependence of the tree contributions A n ^ K (fi), which are 0(p 2 ), is compensated 
by the collections lj %j k of all chiral logarithms, which read 

Ml 



x 



576tt 2 F 2 

(74 + 22r)x 2 + (253 + 88r + 31r 2 )xy - (42 + 12r 2 )xz] ln(/i 2 /A 2 ) , 

1 



uj K (x,y,r) 



+ ^(3x + 8y) + 2/2 K [x + (r + 2)y] + -fi v [x + Ay + 4(r - l)z\ 
Ml 



192vr 2 F2 
1 



x 



(59 + 37r)x 2 + (52 + 45r + 27r 2 )xy - (16 - lOr + 12r 2 )xz 



ln(/x 2 /A 2 ^ 



+ -^{2x + (5 + r)y} + -ii K [2x + 3{r + l)y} 



12 



/^[10x + (5 + I7r)y + 8(r - l)z) 



(5.34) 



In order to get the full expressions of the parameters a and (3 as functions of r up to and 
including the 0(m 2 ) terms, it is sufficient to solve the system of equations ( [5.33D for x and y. 
This can be done perturbatively: x, y and z in the functions u n , ujk can be replaced by their 
leading order values xq, yo, Zq, 



x = 1 - y Q , y = e(r) 



2 2(r- IV 



(5.35) 



since they multiply terms of order 0(p 2 ). The final expression for a as a function of r then 
reads 



a(r) = 1 + 3e*(r) 



4 (Fl 



r - 1 \F2 



1 + 18(2 - r)p! - 6rp 2 + a {2) {r) + 0(p 



(5.36) 



where e*(r) is defined in Eq. ( |5.12|) and the 0(p 2 ) part is given by a^ 2 ^ = «tree + ^loop; with 



a 



(2) 
tree 



6 



a 



(2) 
loop 



r — 1 
6 

r — 1 



(Ax - A,) + A a + A x 



(5.37) 



[ujk(xo, Vo, r) - u„(x Q , y , r)] + uj a (x , y , r) + uj n (x , y , r) . (5.38) 



39 



The various contributions to a and /3, displayed in Eqs. ( |5.36| ) and ( |5.31| ), are shown 
in Figures 8 and 9, respectively. The leading contribution is given by Eq. (|5.8|) , (/^ 
= 1). The O(p) correction consists of terms proportional to F^/F^ — 1 and of the p\ 
and p2 terms, which are enhanced by factors of r. Using the estimate Q5.22 ), the uncertainty 



in + induced by the latter terms is shown as the doubly-hatched area in Fig. 8a. 
Notice that (3^ + (3^ (Fig. 9a) is free from such an uncertainty - cf. the first two terms 
in Eq. ( |5.31| ). Considering the 0(p 2 ) components and (3^ 2 \ two facts should be kept in 
mind. First, there is at present no direct way to determine the constants of £(2,2) an d £(o,4) 
entering the corresponding tree contributions a[f ee (/j,) and Arcc(/-0- Next, the splitting of 
and (3^ 2 ' into the (//-dependent) tree and loop parts is a matter of convention. At present, the 
only available information on and (3^ that is unambiguous concerns the variation of the 
loop contributions (//) and (3^l p (/j) with the renormalisation scale /x (see Figs. 8b and 
9b). We use this variation between \x = M n and \l = 1 GeV as an estimate of the size of the 
contributions of and (3^ . The hatched region of Figs 8a and 9a then shows the total a 
and (3 respectively, with uncertainties arising from the p's and from the 0(p 2 ) contributions 
added. (It should be noted that the final a(r) and (3{r) displayed in Figs, la and lb of Ref. 
26] have been obtained in a slightly different manner: the contributions a\ oop {[i) and /3 loc , p (/i) 



have been simply added to + and (3^ + (3^ respectively. While in the case of a the two 



procedures practically coincide, the curve f3{r) shown in Fig. 1 of Ref. [26] should be considered 
as slightly overestimated.) The final uncertainties in a(r) and j3(r) could be reduced when an 
independent estimate of the constants from £(2,2) and £(0,4) becomes available. 

5.4 The SU{2) L x SU(2) R GOR Ratio 

The above expression for the parameters a and (3 as functions of the quark mass ratio r = m s /m 
is characteristic of the SU(3)l x SU(3)r and it is subject to the corresponding uncertainty. 
One may ask whether there is a direct relation between a and (3 and the condensate defined in 
the SU(2)l x S77(2)r symmetry limit, which would not be based on an expansion in m s of a 
quantity such as M\. The existence of such a relationship can be indeed be inferred already 
from the standard SU(2)^ x SU(2)r one-loop analysis of Gasser and Leutwyler |2B], ffl. In this 



case, the parameters a and (3 are given by Eqs. (|3.29|) and ( 3.3C ) in terms of the constants 
£3 and I4. It turns out that £3 also measures the deviation from the S77(2)l x SU(2)n GOR 
relation [HI: 



2mB M, 



2 



Ml 32vr 2 FJ 

7T 

where 



1 + t^-^4 + 0(M w 4 ) , (5.39) 



B = - lim {Q\uu\Q)/F;. (5.40) 

m s fixed 
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Eliminating the constants £3 and £4, one obtains 

2mB a- 3 Ml A , , 

This formula clearly shows that the ratio on the left hand side should be indeed considerably 
less than 1, provided a is as large as indicated by the central value of existing data. However, 
Eq. ( |5.41| ) is based on the standard and there is no need for it to hold if either a or B 
differ from 1 by more than an amount of order 0(M 2 /A 2 H ). It is convenient to rewrite Eq. 
(|5.41| ) in terms of the GOR ratio 

2m 2mBF 2 
xgor = -j^ Jmjn\uu\n) = -j^r , (5.42) 

m s fixed 

where F = lim mti)md ^ F n . Using the standard xPT relation 

F 2 M 2 

Eq. ( |5.41|) can be reexpressed as 

r -2 a + 2P 3M ' (^44) 
xgor - 2 - — g— - . (5.44) 

Within the standard expansion, Eqs. ( |5.41|) and (|5.44|) are strictly equivalent. It turns out, 
however, that Eq. (|5.44 ) is of a more general validity: It remains true at the leading order of 
G^PT, i.e. irrespective of how much a does deviate from 1. Here, we are going to discuss the 
O(p) and 0(p 2 ) corrections to Eq. (|5.44j ). 

The SU(2) L x S77(2) R condensate ( |5.40|) can be obtained from Appendix A as follows 

F 2 - F 2 M 2 

-^B= lim -2— £ . (5.45) 
F 2 ^ 2mF 2 1 ; 

m s fixed 

The result is given in the form of an expansion in powers of m s , 

F 2 ( 1 

—2 B = B + 2m s Zl + m 2 (p i + ~p 5 + 6p 7 
t< v Z 



+ m] (if + if + 4if s + 2if s 

- 2p K + 2m s A + Qm s Z^ 

- ^ (B + 10m s Z 5 - 8m s Z p ) , (5.46) 
where Mr and M v denote the kaon and eta masses in the SU{2)\ J x SU(2)n chiral limit, and 

m „ = hm = _j In . (5.47) 

m s fixed 
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The formula ( |5.46| ) collects Zweig-rule violating contributions to the SU(2) L x 5?7(2) R conden- 



sate B which are induced by the non- vanishing strange quark mass. It includes, in particular, K 
and 1] loops, which, being proportional to a power of m s (or of r), may be rather large and sen- 
sitive to the subtraction scale /i. This does not mean, however, that the difference F 2 B — FqBq 
of the two-massless flavour and three-massless flavour condensates should be expected to be 
large, but rather that the chiral logarithms alone can hardly be used here in order to estimate 
the size of this difference. 

One may now proceed as follows. First, one subtracts from F 2 M 2 a r given by Eq. (|2.18|) 
the expansion of 4F 2 M 2 displayed in Appendix A. This eliminates fn 2 A . Then, one replaces 
FqBq by F 2 B, using Eq. (|5.41|) . One finds that the dangerous Zweig-rule violating terms 
enhanced by higher powers of r cancel. Next, one eliminates m(£ + 2£) in favour of j3 r , using 
Eq. ( |2.17| ). Finally, one replaces a r and f3 r by a and /3, introducing the unitarity logarithms 
displayed in Eqs. (|3.23|) and ( |3.24|) . As usual (c/. the discussion in Section 11 of Ref. ||), the 



remainders of K and r\ loops are absorbed into the 0(p 4 ) tree-level constants. In this way, one 
arrives at the desired relation 

xgor = 2 - Q + + 15pi -p 2 + Ax(/i) + ZgoVaO • ( 5 - 48 ) 

Here, Ax(fj) = 0{m 2 ) collects all (redefined) tree contributions arising from £(2,2) and £(o,4)- 
The 0(p 2 ) pion-loop contribution reads 

x loop - Ml ( 3Uy I 22^ ln^ 

Xgor ~ 192tt^ I 3 + 22 V ln A 2 

-fi n (l7y + lly 2 ) - (3 + Uy + lly 2 ) , (5.49) 

where the parameter y may be replaced by its leading order expression in the 577 (2) l x 5^7 (2) r 
expansion: 

a + 2(3 . . 

V = ~2Z = 1 - XGOflllead = ~ 1 • (5.50) 



Eq. QQ8D gives the SU{2) L x SU{2) R GOR ratio (^42|) as a function of (a + (3)/3. It is seen 



that for a and (5 close to 1, one recovers the standard formula (|5.44|) up to and including 
0(p 2 ) accuracy. Notice that, using the estimate ( |5.22[ ), the O(p) contribution 15pi - p 2 to xqor 
becomes negligible in the whole range of r. Fig. 10a shows the GOR ratio (|5.48|) with the 0(p 2 ) 



uncertainty obtained from x 1 q°q R {^) as described at the end of Subsection 5.3. 



6 Summary and Conclusions 



In the present article, we have obtained the explicit expression of the tttt scattering ampli- 
tude A(s|t, u) to two-loop accuracy in the chiral expansion. This constitutes our main result, 
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and it is contained in Eqs. ( |3.46| ), ( |3.47| ) and ( |3.49| ). At this order, the Tin amplitude is 
entirely determined by the general requirements of analyticity, crossing symmetry, unitarity 
and the Goldstone-boson nature of the pion, up to the six independent parameters a, (3, Aj, 
i = 1,2, 3, 4, which are not fixed by chiral symmetry. 

For the four constants Aj, we have established a set of sum rules which are evaluated 
using available information on tttt interaction at medium and high energies (E > 600 MeV). 
Details about the derivation of these sum rules and the corresponding numerical analysis will 
be published separately. The resulting values of the A,'s are given in Eqs. ( [4.9|) and ( }4.10|) . For 
Ai and A2, they are compatible with previous determinations, but with a much smaller error 
bar in the case of Ai. The constants A3 and A4 are genuine 0(p 6 ) contributions. To the best of 
our knowledge, they have not been discussed before. 

The most important improvement when going from one loop to two loops concerns uni- 
tarity: The latter is only satisfied in the perturbative sense by the amplitude A(s|£, u). A study 
of the Argand plots for the lowest partial waves {cf. the partial wave f$ in Figure 1) shows 
that violations of the unitarity bounds become important only at energies greater than 450 
MeV. From the size of the two-loop corrections, we also expect that higher orders are under 
control. Typically, for e.g. the S-wave scattering lengths, the one- loop corrections represent a 
25% effect. The two-loop contributions, however, do not modify the one-loop result by more 
than 5%, independently of the value of a. (At this level of accuracy, isospin-breaking effects 
due to rrid 7^ m u and to electromagnetism could start to play a role. We intend to discuss these 
effects elsewhere.) 

As far as the two remaining parameters, a and (3, are concerned, let us stress once more 
the particular importance attached to the first one: a measures the amount of explicit chiral 
symmetry breaking in the irn amplitude due to the non-vanishing of the light quark masses. 
As such, it is intimately correlated to the ratio x = —2m(0\qq\0)/F%M% or equivalently, to the 
quark mass ratio r = m s /m. The commonly accepted picture that spontaneous breakdown 
of chiral symmetry results from a strong condensation of quark-antiquark pairs in the QCD 
vacuum requires that x ~ 1 (or r ~25), and it is only compatible with values of a and (3 close 
to unity, as predicted by standard x?T. A sizeable deviation of a from 1, say e.g. a > 2, even 
if not "a major earthquake" |k3], would nevertheless shake our confidence in the validity of the 
standard lore. These connections between a and the ratios x or r are well established at leading, 
0(p 2 ), and at next to leading, 0(p 3 ), orders, cf. Eqs. (|5.48| ), ( |5.36l ) and (|5.31|) and Figs. 8, 9, 
and 10. Our estimates of 0(p 4 ) corrections suggest that these correlations do not suffer major 
changes even at this order. In the particular case of the standard (r ~ 25), higher order 
corrections to a are under even tighter control: The large-condensate hypothesis can hardly 
be compatible with a > 1.2. Thus, an accurate determination of a and (3 would confirm or 
contradict the standard picture of the chiral structure of the QCD vacuum. Unfortunately, the 
result of a fit to the presently available K14 data |39[] remains inconclusive in this respect. We 
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obtain 

a = 2.16 ±0.86, (3 = 1.074 ± 0.053 . 

Clearly, additional information, which would make such fits more accurate is needed. This may 
be provided by new high statistics Kn experiments, that could be done e.g., at DA$NE, or by 
a precise measurement of the lifetime of 7r + 7r~-atoms at CERN. 

It is, however, interesting to stress that the above fit result leads to values for the thresh- 
old parameters which are in perfect agreement with the results obtained some time ago from 
analyses of data based on the Roy equations. For instance, the fitted value of the 1=0 
S-wave scattering length, a[] = 0.263 ± 0.052, reproduces the result quoted in f24|. Similar 



conclusions hold for the other threshold parameters, see Table 1. It seems thus quite reason- 
able to conclude that in the energy range accessible in decays, our two-loop chiral result 
contains all the relevant information on the irn interaction which is already encoded in the Roy 
equation (for instance, the constraints coming from the production data at higher energies are 
satisfactorily reproduced by our determination of the Aj's). 

The fact that the two-loop expression of the tttt scattering amplitude was obtained, in 
Section 3, in a rather direct and easy way from general S-matrix properties, and without any 
reference to the effective lagrangian C eS or to Feynman diagrams, should not be surprising. Af- 
ter all, the effective lagrangian is merely a technical device which collects the general properties 
of transition amplitudes among Goldstone bosons. All the results obtained with the help of the 
effective lagrangian are, in principle, reproducible using directly Ward Identities, analyticity, 
crossing symmetry and unitarity within a systematic low-energy expansion of QCD correlation 
functions. It just happens that in the case of the 7T7t amplitude, it is technically simpler to 
use the direct S-matrix method than to evaluate all two-loop graphs generated by the effective 
lagrangian. In any case, the results have to be identical in both approaches. We have shown 
that this is indeed the case at the 0(p 4 ) level (see Section 3.2), and since we understand that 
work towards an effective lagrangian computation of ^4(s |t, u) to order 0(p 6 ) in the standard 
framework is in progress |54|, [5f| , we hope that it will soon become possible to check this state- 
ment to two loops by a direct comparison at least in this particular case. Undertaking a similar 
enterprise in the generalised case would be considerably more difficult, due to the proliferation 
of low-energy constants in C eS to that order. For the same reason, it is not clear that it would 
be particularly useful. 
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Appendix A Expansions of Masses and of Decay Constants 

to Order 0(p 4 ) 

From the expression of the effective action at order 0(p 4 ) one computes 



-| = 1 + 2m[£ + (2 + r)i] + 2m 2 5%$ - 4^ - 2fi K 



(A.l) 



and 



|| = 1 + m[(r + l)f + 2(2 + r)£] + 2m 2 5 { ^ - \^ - 3fj, K - \^ , 



where 5p^ and <5j^ j contain the contributions from £(2,2)) c/. Eq. ( |2.12|) 



:(2,2) 



(A.2) 



= A 1 + -A 2 + 2A 3 + -(2 + r 2 )(A 4 + 2B A ) + 5 X - B 2 + 2(2 + r)L>' 



r(2,2) 



= -(1 + r 2 )(A 2 + A 3 + Si) + - (A 2 + 2A 3 - 2S 2 ) 
+ \(2 + r 2 )(A A + 2B i ) + (r + l)(2 + r)D s . 



Similarly, the expansion of the pion and kaon masses read 

2mB + 4m 2 A + 4(2 + r)m 2 Z$ 



F 2 

—M 2 



2m 3 <J + 2m 4 5^ + 4m 2 M 2 A 3 



4r(0,4) 



3M 2 + 20m 2 (A + 2Z^ 



2fx K M 2 + 4(l + r)m 2 (A + 2Zo' 



Ml + 4(1 + 2r)m 2 {A Q + 2Z 5 ) + 8(1 - r)m 2 {A + 2Z P ) 



(A.3) 



(A.4) 



Ek m 2 



= (r + l)mB Q + (r + l) 2 m 2 A + 2(r + 1)(2 + r)m 2 Z^ 

+ (r + l)m 3 5g'| + (r + l)m 4 <'| + (r + l) 2 m 2 M|A 3 
3 



Mr- + 4(r + l)m 2 (A + 2^ 



3/i K + 2(1 + r) 2 m 2 (A + 2Z £ 
1 

6^ 



5M£ + 4(r + 1)(1 + 2r)m 2 (A + 2Z$) - 4(1 - r 2 )m 2 (A + 2Z, 



(A.5) 
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The 0(p 3 ) tree-level contributions 8^ 3 p, P = ir, K are given as follows 
<4J = ^i + ^2 + (10 + 4r + r 2 )p 4 + i(2 + r 2 )p 5 + 6(2 + r) 2 p 7 , 

(A.6) 

= ^ + r + r 2 ) Pl + 1(1 - r + r 2 )p 2 + 3(2 + 2r + r 2 )p 4 + 1(2 + r 2 )p 5 + 6(2 + r) 2 p 7 . 

The 0(p 4 ) tree level contributions from £(°' 4 ) contained in and in 5^'^- are not displayed 
explicitly. 
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Appendix B The Polynomials £^(s) 



We give here the list of polynomials ^'(s) which define the one-loop partial- wave pro- 
jections in Eqs. ( |3.34| ) and (|3.36| ). The center-of-mass momentum q of the pions is given (in 
units of M n ) by 



s-4M 2 



(B.l) 



+ 



1447T 2 

1 

72^2 
311 



35a 2 + 80a/3 + 134/3 2 J + 10 (A x + 2A 2 
(60a + 209/?)/? + 16(2Ai + 3A 2 )| q 2 



1 



1927T 2 

1 

11527T 2 

1 

288tt 2 



1 



7 



5a + 16(3 + 24f3q 2 
1 



2 „4 



-5a 2 + A(3 Z 



+ 



0V 



12tt 2 



(B.2) 



+ 
+ 



+ 
+ 



1 



288tt 2 
1 



1447T 2 



f 265 



12167T 2 

1 

576vr 2 
1 



31a 2 - 122a/3 + 220/3 2 J + 4 (A a + 2A 2 

(-69a + 268/?)/? + 8(Ai + 3A 2 )| q 2 

1 6 1 
/5 2 + y(A 1 + 4A 2 )]g 4 



1447T 2 

11 



9a 2 - 42a/3 + 60/? J 
-9a + 37/3) /3 g 2 



2 „4 



727T 2 



(B-3) 
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d 2) ( s ) 
d 3) ( 5 ) 
d 4) ( 5 ) 



288tt 2 
1 

288tt 2 



a -4(3- 6(3 q 2 
-3a 2 + 2a(3 - 12(3 2 



24tt 2 



= 



^°\s) = ^ {5a 2 - 80aP + 10P 2 } 

J _JL_ [ 55a _ 68/3] /3 - | (Ax - A 2 )| g 2 



+ 



d 2) ( 5 ) 
d 3) ( 5 ) 
e! 4) (.) 



14327T 2 
108tt 2 ' 3 



2 + ? (Ai - A 2 ) f g 4 



1 



288tt 2 
1 



-5a + 7(3) (3+ {5a -3(3)(3q 2 -^ (3 2 q 4 



1447T 2 



72tt 2 



72tt 2 
1 

864tt 2 
5 



144vr 2 



-5a 2 + 10a/3 + 28/3 2 
a 2 + 4a/3 - 2/3 2 



+ 



2 „2 



247T 2 



(B-4) 
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Appendix C The Loop Integrals K n (s) 



In this Appendix, we discuss some properties of the functions K n (s), cf. Eqs. ( |3.42 ), 
( 3.37 ) and ( |3.49|) . Let us start with the function J(s) which denotes the usual <i-dimensional 
one-loop integral 

' (CI) 



(2nY (M*-?)[M*-(q-p)*] 



subtracted at p 2 = 0. J(s) is finite for d = 4 and may be expressed in the standard parametrical 
form [@, | 

J(a) - ' 



16tc 2 Jo 



d\ In 



(C.2) 



For s < 4M 2 , the latter formula can be integrated by parts and, after introducing a new 
integration variable 

Ml 



x 



A(l-A) ' 

it can be transformed into the dispersive representation ( 3.14 ). For real s, one has 



2 + <7(lnf^+i7r) 

' AMl-s\ 2 



if s > 4Mf 



l^J(s) = { 2-2(^) 5 arctg( n ^) 5 , if < a < AM* 



2 + alnf-i, 



if s < 0, 



(C.3) 



(C.4) 



where 



a 



1 - 



4M 



2\ V2 



(C.5) 



The function J(s) was known at times as the Chew-Mandelstam function p6 . 

In order to simplify the discussion of the two-loop integrals, it is convenient to introduce 
the notation 

F{s) = 16tt 2 J(s) - 2 . (C.6) 

The function F(s) vanishes at threshold; more precisely, for s — 4M 2 —> + , one has the 
expansion 



F(s) = ina-2a 2 " ~ 4 



—a <t + . . . 

3 5 



whereas for s — > 0, 



(C.7) 



(C.8) 



The dispersion integrals ( p.42| ) can now be directly determined by algebraic manipulations with 
the functions F(s), without calculating a single integral. It is sufficient to construct functions 
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that are analytic except for a branch-cut singularity on the positive real half-axis starting at 
s = AM 2 , with the discontinuities given by the functions k n (s), Eqs. ( 13.371) : 



lmK n (s) = ^[Kn(s + ie) - K n (s - ie)} = k n (s)6(s - AM 2 ) . 



(C.9) 



Such analytic functions are, of course, only defined up to a polynomial. The latter, however, 
can be unambiguously fixed, using the boundary conditions 



lim s^K^s) = , K n (Q) = Q, 

s|^oo 



(CIO) 



which follow from Eqs. ( |3.37| ) and ( |3.42| ). The main ingredients of this construction are the 
analytic functions F(s), F 2 (s), and F 3 (s), with the discontinuities 



ImF(s) 
lmF 2 (s) 

lmF 3 (s) 



2tt 



s-AMl 



L{s) 6{s - AMI) 



s-4M. 



2\ 3 / 2 



(3ttL 2 (s) -ti 3 )9(s-4M- 



7r) ' 



(C.11) 
(C.12) 

(C.13) 



where L(s) is the logarithmic function ( 3.38Q . Since F(4M 2 ) = 0, new analytic functions with 
the desired discontinuities can be obtained upon dividing by (s — 4M 2 ). This leads to 



Im 



F(s) 



TT 



s - 4M ' Js (s - 4M2) 



0{a - AMI) 



(C.14) 



Im— ±L = 2kL{s)6{s-AMI) 



(C.15) 



This procedure can be extended further: Subtracting from the function sF 2 (s)/(s — AM 2 ) its 
threshold value 

sF 2 (s) 



lim 



4M2 S - AM 2 



-7T 



(C.16) 



one obtains from Eq. (C.15) 



Im 



1 



s-M 2 



sF 2 (s) 
s-M 2 



2ttL(s) 
s-M 2 



(s - AM 2 ) . 



(C.17) 



The functions k (s), ki(s) and k 2 (s) can be recognised on the right-hand sides of Eqs. (C.ll), 
(C.15) and (C.12) respectively. Similarly, Eqs. (C.14) and (C.17) reproduce the first two 
terms in the expression of k±(s). The remaining terms, quadratic in L(s), can be inferred from 
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Eq. (C.13): The function F 3 (s) + vr 2 (s - AM 2 )F(s)/s exhibits a double zero at s = AM 2 . 
Consequently, one can write 



Im 



> - 4M 2 ) 2 



7T 



F 3 (s) + — (s-AM 2 )F(s) 
s 



3nL 2 (s) 



- AMI) 



9{s - AMI) ■ 



(C.18) 



Finally, subtracting from the analytic function on the left-hand side of the last equation its 
threshold value, one obtains 

i r mi 



Im 



s - AMI I (s - 4M 2 ) 2 



sF 6 (s) +^(s-AM:)F(s) 



11 



_3ttMI_ 
s-AMl 



L 2 (s) 



■■ 6(s - AMI 



(C.19) 



's{s - AMI) 

All the elements needed to write down the functions K n (s) are now collected. Eqs.(C.ll), 
(C.15), (C.12), (C.18) imply, in succession, 

1 



Ms) 

Ki(s) 



16vr 2 
1 s 



F{s) + Kq(s) 

F 2 (s) + k 1 (s) 



Vo^s-AMT 



16?r 2 
1 



F 2 (s) + k 2 (s) 
sM 2 [ 



(C.20) 



q _ AM 2 

F 3 (s) + n 2 S * F{s) 



16tt 2 (s - AM 2 ) 2 l 
and from Eqs. (C.14), (C.17) and (C.19) one reconstructs K^{s): 
Ml F(s) 



KAs) 



16tt 2 s - AM 2 



+ 
+ 



M 2 



1 



32tt 2 s - AM 2 
Ml 1 



s -AM 2 
Mis 



F 2 (s) + n< 



(C.21) 



487r 2 s-4M 2 {(s-AM 2 y 



F 3 {s)+7c' 



,s-AM, 



-F(s) 



7T 2 \ + K 4 (s) . 



In Eqs. (C.20) and (C.21), « n (s) denote sofar arbitrary polynomials. They may be fixed using 
the boundary conditions (C.10). Since F(s) is bounded by Ins for \s\ — > oo, the asymptotic 
condition (C.10) implies that all K n (s) must be constant. Their values are then determined by 
the conditions K n (0) = 0: 



Kq 



in' 



K\ = 0, K 2 



1 



1 



4tt 2 



132 32tt 2 



(C.22) 



Eqs. (C.20) and (C.21) now coincide with the formula (|3.49 ) given in Section 3 without proof. 
It is rather interesting that all two-loop contributions to the irn scattering amplitude can be 
obtained by algebraic manipulations with the Chew-Mandelstam function, without having to 
evaluate a single integral. 
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Appendix D Threshold Parameters 



In this Appendix, we display the explicit two-loop expressions of the threshold parameters 
aj and b\. They are defined from the low-g 2 expansions of the real parts of the corresponding 
partial waves //(s) 

Refi(s) = q 21 {a] + b]q 2 + ■■■} , (D.l) 

with q 2 = - AMI) /AMI. The partial waves are defined by taking the corresponding projec- 
tions O of the two-loop amplitude ( gg§ ), (|3^ ). 

1 M 2 5 M 4 1 /If 4 

^(5a + 16/3) + --HA! + 2A 2 ) + — — --£(5a + 16/?) 2 



96vr Fy 87r F 4 y 4608vr 3 F 4 

1 M 6 5 M 6 

7T 7T 



1 /? 1M 2 . oA 1 M 2 {3 2 5 5 



2^ 

4vr FJ vr Fi v 7 3vr 3 Fi v 3 96 r 256 



7T 7T 7T 

5 M i „«, 81 » i 5 K ,13 37 .. .„ 
~ 33T77faS7f (161+ 24 T)a + 92T6^7f ( T-lT T)a ' 3 

+ 5 ^ 85 1^ 4793 _ 823 (D . 3) 

20736tt 5 F? v A 1 H 41472tt 5 Fj? 2 3 JH K 1 



5vr F6 v y 36tt 3 F6 1V 72 225^ 36vr 3 F^ zv 36 225 r 

11 Ml 41 9x , 1 M 2 121 9N „„ 

+ ;r^r 7T a 7i 2 )a 2 6 

7464960tt 5 F6 v 4 ; 138240vr 5 F^ 18 ' H 

1 M 2 46 1 M 2 1679 

U - '39 - -tt>/3 2 + ( — - 3297T 2 )/? 3 (D.4) 



248832vr 5 F5 v 5 ' 9331200tt 5 F2 



52 



1 /W 4 1 M 4 5 1 /If 4 5 

+ 2^7f A *> + IS^Tf A '<I2« " » + sT^Tf A %« - « 

5 M 4 7T 2 3 5 M 4 7T 2 2 

+ 497664vr 5 F6 ( 6 j " + 331776tt 5 F% { 3 j " 1 

5 M 4 ,374 in 2N n2 5 M 4 73 r 2 . - /T ^ rX 



497664tt 5 F6 v 3 746496tt 5 FJ? v 2 



* = -^ (Al - A2) + 4320^^Y /3 a -T^"24^ 

+ ^ (As " A4) + 48^^ Al( 2 ~ 15^ + W A2( 2 ~ 
1 ,107 2 , 3 1 M l /H3 2 , 2a 



746496tt 5 F6 v 8 ' 27648vr 5 v 36 

1 Ml ,851 47 , 2 1 ,4601 589 2 3 
+ 41472^ ( T" ~ T ^ " 1244160^^ ( ^ 3^ )(3 (D ' 6) 



1 1 23 1 211 

+ «(A 3 -A 4 )H — -Ai( a-/3)H — — A 2 (a 13) 

35ttF6 v 3 4; 294tt 3 F6 u 240 p; 3528tt 3 F6 2V 10 1 ' 

1 .163 2 o 1 311 2 , 2 

+ ( 7T )a H f 7T W /3 

20321280tt 5 F6 1 12 ; 8 1285 12^^ 1 10 ; ^ 

1 ,2537 „ 9 x ^ 1 ,8011 13 

(— - 17tt 2 )^ 2 + — — - — (— - - — tt 2 )/? 3 (D.7) 



2903040tt 5 F6 v 14 ' ^ 3386880tt 5 F6 v 36 5 



1 M T 6 , 1 MJ? , , 

^ A 3 + T7-^(Ai + 2A 2 )(a-4/?) 



2tt F6 ° 48tt 3 F% 

1 M? ,29 9x , 1 Ml , 23 2s 9^ 

+ i»^(y-^> -9^^f (37 -y^ 

+ 1536^^ (1? " 18" ^ " 384^ HJ ( T " Y )P (D ' 8) 
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6 » = ^ir : + h^ +3X ^^M {112p2 - 8ap - 7a2) 



2tt Fy ° *' 432vr 3 F6 1V 4 ^ 108vr 3 F% 

1 "l(n-?,V-^^( 9 i-^V^ 



165888vr 5 Fy 6 ' 36864vr 5 i 7 ; 6 v 9 
+ 82944^7* (1487 " ~2~ M " 82944^7* ( ^ 3^ )/5 (D ' 9) 



2 30ttF4 v 27 7200^^8 6 P 3 P; 

1 M 2 , x 1 M 2 151 . 1 M 2 41 

£ ( A 3 + A 4 ) - onAnn Ai(— a - 191/5) - -^^-^ A 2 (— a - 61/5) 



5tt F6 v J 4y 32400tt 3 F6 1V 2 ^ 8100tt 3 i^ 6 zv 2 

1 M -- 1223 - 7T y + 7 M ( lH_, V g 



9331200tt 5 F6 v 16 y 1382400tt 5 i^ 6 v It 

1 K% 11n 257 2 2 1 M 2 3289 , 3 

+ 518400^ Ff (U9 " "24 ^ + 18662400^1* ( ~ " 271?r )/5 (D ' 10) 

The contributions coming from various orders of the chiral expansion are easy to distin- 
guish. The lowest order contributions, corresponding to the amplitude (|3.9| ), give terms linear 
in a and (3, proportional to F~ 2 . At this order, only the S- and P-wave scattering lengths and 
the S-wave slope parameters are non-vanishing. The one-loop precision brings in two types of 
contributions, both proportional to F~ A : Those which are linear in Ai and A 2 come from 0(p 4 ) 
and 0(p 5 ) tree graphs, and those which are quadratic in a, (3 correspond to genuine one-loop 
graphs. Finally, the contributions at order two-loop, proportional to F~ 6 , come from 0(p 6 ) 
and 0(p 7 ) tree-graphs (terms linear in A3 and A4), from one- loop graphs with a Ai or A 2 vertex, 
whereas the genuine two-loop graphs generate the contributions which are cubic in a, (3. 



54 



References 



[1] G. 't Hooft, in: Recent developments in gauge theories, eds. G. 't Hooft et al. (Plenum, 
NY, 1980); Y. Frishman, A. Schwimmer, T. Banks and S. Yankielowicz, Nucl. Phys. B177 
(1981) 157; S. Coleman and B. Grossman, Nucl. Phys. B203 (1982) 205. 

[2] C. Vafa and E. Witten, Nucl. Phys. B234 (1984) 173; Comm. Math. Phys. 95 (1984) 257. 

[3] S. Weinberg, Physica 96A (1979) 327. 

[4] J. Gasser and H. Leutwyler, Phys. Lett. B125 (1983) 321; Ann. Phys. 158 (1984) 142. 

[5] J. Gasser and H. Leutwyler, Nucl. Phys. B250 (1985) 465. 

[6] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; 124 (1961) 246. 

[7] M. Gell-Mann, R. J. Oakes and B. Renner, Phys. Rev. 175 (1968) 2195; S. Glashow and 
S. Weinberg, Phys. Rev. Lett. 20 (1968) 224. 

[8] N. H. Fuchs, H. Sazdjian and J. Stern, Phys. Lett. B269 (1991) 183. 

[9] J. Stern, H. Sazdjian and N. H. Fuchs, Phys. Rev. D47 (1993) 3814. 

[10] B. Adeva et al, "Lifetime measurement of 7r + 7r~ atoms to test low-energy QCD predic- 
tions," preprint CERN/SPSLC 95-1 (Dec. 1994). 

[11] M. Baillargeon and P. J. Franzini, "Accuracies of Parameters at DA$NE," contribution 
to the Second DA$NE Physics Handbook, L. Maiani, N. Paver and G. Pancheri, eds., May 
1995. 

[12] P. W. Anderson, Basic notions in condensed matter physics" , Frontiers in Physics, Lecture 
Notes Series, The Benjamin/Cummings Publishing Company, Inc., 1984. 

[13] T. Banks and A. Casher, Nucl. Phys. B168 (1980) 103; E. Floratos and J. Stern, Phys. 
Lett. B119 (1982) 419; D.I. Dyakonov and V.Yu. Petrov, Phys. Lett. B147 (1984) 351; 
Nucl. Phys. B245 (1984) 259; A.V. Smilga and J. Stern, Phys. Lett. B318 (1993) 531. 

[14] H. Leutwyler and A. V. Smilga, Phys. Rev. D46 (1992) 5607. 

[15] J. Stern, in preparation. 

[16] R. Gupta in LATTICE 94, Nucl. Phys. B (Proc. Supp.) 42 (1995) 85. 

[17] M. Knecht and J. Stern, "Generalized Chiral Perturbation Theory", contribution to the 
Second DA$NE Physics Handbook, L. Maiani, N. Paver and G. Pancheri, eds., May 1995; 
preprint IPNO/TH 94-53, |Eep^pTi/94Tl253. 



55 



[18] S. Weinberg, in: A Festschrift for I. I. Rabi, ed. L. Motz (New York Academy of Sciences, 
New York, 1977), p. 185. 



[19 
[20 
[21 
[22 

[23 
[24 
[25 
[26 

[27 
[28 
[29 

[30 

[31 
[32 

[33 
[34 
[35 
[36 



M. Knecht, B. Moussallam and J. Stern, Nucl. Phys. B429 (1994) 125. 
N. H. Fuchs, H. Sazdjian and J. Stern, Phys. Lett. B238 (1990) 380. 
R. Koch and E. Pietarinen, Nucl. Phys. A336 (1080) 331. 

T. E. O. Ericson et ai, "The irNN coupling from high precision np charge exchange at 
162-MeV," preprint CERN-TH-95-50. 

J. Gasser and H. Leutwyler, Phys. Lett. B125 (1983) 325. 

M. M. Nagels et ai, Nucl. Phys. B147 (1979) 189. 

S. Weinberg, Phys. Rev. Lett. 17 (1966) 616. 

M. Knecht, B. Moussallam and J. Stern, "The tttt Amplitude in Generalized Chiral Per- 
turbation Theory", contribution to the Second DA$NE Physics Handbook, L. Maiani, N. 



Paver and G. Pancheri, eds., May 1995; preprint IPNO/TH 94-54, |hep-ph/ 94 11259 
H. Leutwyler, Ann. Phys. 235 (1994) 165. 
M. Knecht, in preparation. 

H. W. Fearing and S. Scherer, "Extension of the Chiral Perturbation Theory Meson La- 
grangian to Order p 6 ", TRIUMF Report TRI-PP-94-68 (Aug. 1994), and |hep-ph/940834e . 

V. Bernard, N. Kaiser and U.-G. Meifiner, Nucl. Phys. B357 (1991) 129; Nucl. Phys. 
B364 (1991) 283. 

J. Bijnens, G. Colangelo and J. Gasser, Nucl. Phys. B427 (1994) 427. 

J. Gasser, "The tctt Scattering amplitude in Chiral Perturbation Theory", contribution to 
the Second DA$NE Physics Handbook, L. Maiani, N. Paver and G. Pancheri, eds., May 
1995, and |hep-ph/941239l 

M. G. Olsson, Phys. Rev. 162 (1967) 1338. 

B. Hyams et al, Nucl. Phys. B64 (1973) 134. 

W. Ochs, vrN Newsletter 3 (1991) 25. 

D. Morgan and M. R. Pennington, "Low Energy Scattering" , contribution to the Second 
DA$NE Physics Handbook, L. Maiani, N. Paver and G. Pancheri, eds., May 1995. 



[37] M. R. Pennington and J. Portoles, Phys. Lett. B344 (1995) 399. 



56 



[38] J. Gasser and U.-G. MeiBner, Phys. Lett. B258 (1991) 219. 

[39] L. Rosselet et al, Phys. Rev. D15 (1977) 574. 

[40] D. Morgan and G. Shaw, Nucl. Phys. BIO (1969) 1387. 

[41] J.-L. Basdevant, CD. Froggatt and J. L. Petersen, Phys. Lett. B41 (1972) 173, 178. 
[42] J. L. Petersen, CERN Yellow Report 77-04 (1977). 

[43] A. Manohar and H. Georgi, Nucl. Phys. B234 (1984) 189; H. Georgi and L. Randall, Nucl. 
Phys. B276 (1986) 241; H. Georgi, Phys. Lett. B298 (1993) 187. 

[44] C. A. Dominguez and E. de Rafael, Ann. Phys. 174 (1987) 372. 

[45] J. Bijnens, J. Prades and E. de Rafael, Phys. Lett. B348 (1995) 226. 

[46] J. Stern, N.H. Fuchs and M. Knecht, "Light quark Masses from Exclusive r Decays: An 
Experimental Proposal", in Proceedings of the Third workshop on the Tau-Charm Factory, 
Marbella, Spain, 1-5 June 1993", J. Kirkby and R. Kirkby, eds. (Editions Frontieres, 1994). 

[47] M. Jamin and M. Miinz, Z. Phys. C66 (1995) 633. 

[48] K. G. Chetyrkin, C. A. Dominguez, D. Pirjol and K. Schilcher, Phys. Rev. D51 (1995) 
5090. 



[49] S. Narison, Montpellier preprint PM-95-06 (April 1995) and |hep-ph 9504333 . 
[50] H. Leutwyler and M. Roos, Z. Phys. C 25 (1984) 91. 

[51] K. Maltman and D. Kochtan, Mod. Phys. Lett. A5 (1990) 2457; J. F. Donoghue, B. R. 
Holstein and D. Wyler, Phys. Rev. D47 (1993) 2089; J. Bijnens, Phys. Lett. B306 (1993) 
343; R. Urech, Nucl. Phys. B433 (1995) 234; H. Neufeld and H. Rupertsberger, Univ. of 
Wien preprint, UWThPh-1995-18 (June 1995). 

[52] J. Gasser and H. Leutwyler, Phys. Rep. 87 (1982) 77. 

[53] H. Leutwyler, "Foundations and Scope of Chiral Perturbation Theory," Berne preprint 
BUTP-94-18, |hep-ph / 9409423| ; to be published in "Chiral Dynamics: Theory and Exper- 



iment," A.M. Bernstein and B.R. Holstein, eds. (Springer- Verlag, 1995) Lecture Notes in 
Physics, 452. 

[54] G. Colangelo, Phys. Lett. B350 (1995) 85. 

[55] J. Bijnens et al, work in progress. 

[56] G. F. Chew and S. Mandelstam, Phys. Rev. 119 (1960) 467. 



57 



[57] P. Estabrooks and A.D. Martin, Nucl. Phys. B79 (1974) 301. 
[58] W. Hoogland et ai, Nucl. Phys. B126 (1977) 109. 



58 



Figure Captions 



1- Argand diagram of the amplitude /° (i.e. y 1 — 4M 2 / s Im/g versus v/l — 4M 2 / s Re/° ) . 
The result at one loop is represented by the dashed line and the result at two loops by 
the solid line. Both cases correspond to the same parameter choice a = 2, (3 = 1.08 and 
\ as in (4.9) and (4.10). 

2- Phase shifts computed at increasing chiral order (see formulas (4. 11)- (4. 13) and the 
discussion following them ). The tree, one- loop and two-loop orders are represented by 
the dotted, dashed and dash-dotted lines respectively. Shown for comparison are the data 
of Ref. H and [§7]] for 5° and Ref. [|8| for <Sg and b\. For 6{ and 6% the Breit-Wigner 



contributions of the p and the f 2 resonances are shown as crosses. 

3- Illustration of the sensitivity of the phase shifts to the values of the parameters a and 
(3. The dotted line is the result of the standard xPT (a = 1.04, (3 = 1.08), the dashed 
line corresponds to a = 2, (3 = 1.08 and the solid line to a — 3, (3 — 1.12. 

4- The phase shift difference 5® — 5\ in the range of energies accessible in the decays 
is shown: a) for increasing chiral orders as in Fig. 2 and, b) for several values of a and 
(3, as in Fig. 3. 

5- Influence of the uncertainties affecting the evaluation of Ai, A4 on the prediction for 
<5q — 5{ . The band delimited by solid lines corresponds to varying the Aj's inside the 
ranges indicated in (4.9) and (4.10) with a, [3 fixed to the standard values, whereas 
the band delimited by dashed lines corresponds to a = 2, (3 = 1.08. 

6- Plot of the correlation between a and (3 as inferred from the Morgan-Shaw universal 
band given by Eq. (^4.14j) . The hatched region covers the set of points consistent with Eqs. 
( |5.36| ) and ( |5.31| ). If the 0(p 2 ) terms are ignored, this set shrinks to the doubly-hatched 



sub-region shown. See the discussion in Section 5. 

7- Predictions for the scattering lengths as functions of a. The inner bands are obtained 
by varying the Aj's inside their error bars with (3 fixed at the center of the Morgan-Shaw 
band. Including the variations in (3 allowed by the latter gives the outer bands. 

8- In a), we show the parameter a(r) at increasing chiral order (see Eqs. (5.36) and the 
following discussion). The doubly-hatched area represents uncertainties in the contribu- 
tions from pi, p2, while the hatched band includes the estimated uncertainty of a^ 2 \ as 

(2) 

described in the text. In b) we show the contribution a l0 Q P (r) for different values of \i. 
Each band corresponds to variations of the scale A in the range M v < A < 1 GeV. 
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9- In a), we show the parameter /3(r) at increasing chiral order (see Eq. ( |5.31 ) and the 



discussion following Eq. (|5.36|) ). The hatched area represents the estimated uncertainty 



in contribution from (3^. In b) we show the contribution (3^l p (r) for different values of 
fi. Each band corresponds to variations of the scale A in the range M v < A < 1 GeV. 

10- In a) we show the Gell-Mann-Oakes-Renner ratio xgor as a function of the linear 
combination (a + 2/3) /3 (see Eq. (5.48)). The hatched area represents the expected 
range of 0(p 2 ) contributions. The vertical dotted lines show the range in (a + 2/3) /3 
corresponding to the fit values Eq. (4.15). In b) we show the contribution x X q°q R for 
different values of /a. Each band corresponds to variations of the scale A in the range 
M v < A < 1 GeV. 
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Argand diagram 
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